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ABSTRACT 


To  determine  and  evaluate  significant  characteristics  of  runway 
roughness,  40  runway  elevation  profiles  from  14  runways  were  analyzed. 
Grade  effects  were  removed  by  high-pass  filtering  at  a  maximum  cutoff 
wavelength  of  500  ft  and  the  resulting  profiles  were  shown  to  be  both 
nonstationary  and  nongaussian.  Since  nonstationarity  was  caused  by  the 
long  wavelength  undulations,  the  profiles  were  again  filtered  at  a  maximum 
cutoff  wavelength  of  100  ft.  These  short  wavelength  profiles  are  shown  to 
be  sufficiently  stationary  for  a  compound  characterization.  The  short 
wavelength  components  are  modeled  in  the  frequency  domain  while  the 
long  wavelength  components  are  modeled  by  the  distribution  of  relative 
maxima  and  minima  per  1000  ft  and  the  joint  distribution  of  wavelength 
and  amplitude  of  the  maxima  and  minima,  A  method  is  presented  for 
generating  a  simulated  runway  which  has  the  characteristics  of  the  real 
profiles  used  in  the  analysis. 


TABLE  OF  CONTENTS 


SECTION  Page 

1  INTRODUCTION .  1 

2  DATA  PROCESSING  . . 6 

2.  1  Description  of  Raw  Data  . . .  6 

2,  2  Filtering  Methods  . .  6 

2.  3  PSD  Computation . .  8 

2.4  Long  Wavelength  Description  . .  11 

2.  5  Other  Descriptors . 12 

3  DATA  ANALYSIS  . 14 

3.  1  Maximum  Wavelength  of  500  Feet  .......  15 

3.  1.  1  Amplitude  Probability  Densities  ....  16 

3.  1.  2  Fower  Spectral  Densities .  23 

3.  1.  3  Stationarity  . .  32 

3.  2  Maximum  Wavelength  of  100  Feet  .......  40 

3.  2.  1  Amplitude  Probability  Densities  ....  42 

3.  2.  2  Power  Spectral  Densities .  46 

3.  2.  3  Stationarity  . .  46 

3.  3  Discrete  Long  Wavelength  Description  .....  48 

4  SIMULATED  RUNWAYS  .  54 

4.  1  Sample  Result . .  55 

4.  2  Other  Uses  of  Simulation  Program .  77 

5  SUMMARY  AND  CONCLUSIONS  . .  78 

APPENDIX 

I  DIGITAL  FILTERING  CONCEPTS . 80 

II  POWER  SPECTRAL  DENSITIES  FOR  TWO 

METHODS  OF  COMPUTATION .  86 

III  POWER  SPECTRAL  DENSITIES .  119 

IV  RMS  HISTORIES  FOR  ALL  LINES  OF  DATA 

FILTERED  AT  500  VFET  . 160 

V  CUMULATIVE  DISTRIBUTIONS  OF  AMPLITUDES  FOR 

ALL  LINES  OF  DATA  FILTERED  AT  100  FEET  ...  168 

VI  RMS  HISTORIES  FOR  ALL  LINES  OF  DATA  FILTERED 

AT  100  FEET .  183 

REFERENCES .  191 


v 


LIST  OF  ILLUSTRATIONS 

FIGURE  Page 

1  Response  of  500  ft  Cutoff  Filter  . .  9 

2  Response  of  100  ft  Cutoff  Filter  ...........  9 

3  Cumulative  Distributions  of  Amplitudes  Filtered 
at  500  ft  for  Travis,  Centerline,  and  George, 

Centerline  . . 18 

4  Cumulative  Distributions  of  Amplitudes  Filtered 

at  500  ft  for  McGuire,  Centerline,  and  Travis,  75L  .  .  19 

5  Cumulative  Distributions  of  Amplitudes  Filtered 

at  500  ft  for  Bakalar,  50L,  and  Torrejon,  12L  ....  20 

6  Histogram  of  RMS  Values  With  Profile 

Filtered  at  500  ft . 24 

7  Summary  Plot  of  Normalized  Spectral  Values .  28 

8  Histogram  of  PSD  Mid  Range  Slopes  . .  30 

9  Power  Spectral  Density  of  Unfiltered  Profile 

from  Altus,  Centerline  .  33 

10  Power  Spectral  Density  of  Unfiltered  Profile 

from  Bakalar,  Centerline  .  34 

11  Power  Spectral  Density  of  First  Third  of  McGuire, 

Centerline .  36 

12  Power  Spectral  Density  of  Middle  Third  of  McGuire, 

Centerline  ,  . .  37 

13  Power  Spectral  Density  of  Last  Third  of  McGuire, 

Centerline . 38 

14  RMS  Histories  for  Altus  and  Bakalar  Filtered 

at  500  ft .  41 

15  Histogram  of  RMS  Profiles  Filtered  at 

Maximum  Wavelengths  of  100  ft  . .  ,  45 

16  RMS  Histories  for  Altus  and  Bakalar  Filtered 

at  100  ft  .  47 

17  Cumulative  Distribution  of  Primary  Peaks 

per  1000  ft . 50 

18  Cumulative  Distribution  of  Half  Wavelengths .  52 

19  PSD  from  Simulated  Runway  . .  58 


LIST  OF  ILLUSTRATIONS  (continued) 


FIGURE  Page 

20  Amplitude  Probability  of  Simulated  Runway .  59 

21  Profile  of  Simulated  Runway . .  61 

22  PSD  Altus  75L  (FFT  -  Prewhitened  and  Smoothed)  .  .  87 

23  PSD  Altus  CL  (FFT  -  Prewhitened  and  Smoothed)  ,  .  88 

24  PSD  Altus  75R  (FFT  -  Prewhitened  and  Smoothed)  ,  .  89 

25  PSD  Bakalar  50L  (FFT  -  Prewhitened  and  Smoothed)  .  90 

26  PSD  Bakalar  CL  (FFT  -  Prewhitened  and  Smoothed)  .  .  91 

27  PSD  Bakalar  50R  (FFT  -  Prewhitened  and  Smoothed)  .  92 

28  PSD  Carswell  13L  (FFT  -  Prewhitened  and  Smoothed)  .  93 

29  PSD  Carswell  87L  (FFT  -  Prewhitened  and  Smoothed)  .  94 

30  PSD  Carswell  13R  (FFT  -  Prewhitened  and  Smoothed)  .  95 

31  PSD  Edwards  75L  (FFT  -  Prewhitened  and  Smoothed)  .  96 

32  PSD  Edwards  CL  (FFT  -  Prewhitened  and  Smoothed)  .  97 

33  PSD  Edwards  75R  (FFT  -  Prewhitened  and  Smoothed)  ,  98 

34  PSD  Ellsworth  CL  (FFT  -  rewhitened  and  Smoothed)  ,  99 

35  PSD  Langley  SOL  (FFT  -  Prewhitened  and  Smoothed)  ,  100 

36  PSD  Langley  CL  (FFT  -  Prewhitened  and  Smoothed)  .  .  101 

37  PSD  Langley  50R  (FFT  -  Prewhitened  and  Smoothed)  ,  102 

38  PSD  Altus  7  5L  (Autocorrelation  *  Prewhitened  and 

Smoothed) . , . .  103 

39  PSD  Altus  CL  (Autocorrelation  -  Prewhitened  and 

Smoothed) . 104 

40  I-’SD  Altus  75R  (Autocorrelation  -  Pre'vhitened  and 

Smoothed)  . .  105 

41  PSD  Bakalar  50L  (Autocorrelation  -  Prewhitened  and 

Smoothed)  . .  106 

42  PSD  Bakalar  CL  (Autocorrelation  -  Prewhitened  and 

Smoothed) . 107 

43  PSD  Bakalar  5QR  (Autocorrelation  -  Prewhitened  and 

Smoothed) . 108 


LIST  OF  ILLUSTRATIONS  (continued) 


FIGURE  Page 

44  PSD  Carswell  13L  (Autocorrelation  -  Prewhitened  and 

Smoothed)  . . , .  109 

45  PSD  Carswell  87 L  (Autocorrelation  -  Prewhitened  and 

Smoothed) . * . .  .  .  ,  110 

46  PSD  Carswell  13L  (Autocorrelation  -  Prewhitened  and 

Smoothed) . Ill 

47  PSD  Edwards  75L  (Autocorrelation  -  Prewhitened  and 

Smoothed) . . . . .  1 1 Z 

48  PSD  Edwards  CL  (Autocorrelation  -  Prewhitened  and 

Smoothed)  . .  113 

49  PSD  Edwards  75R  (Autocorrelation  -  Prewhitened  and 

Smoothed) . 114 

50  PSD  Ellsworth  CL  (Autocorrelation  -  Prewhitened  and 

Smoothed) . 115 

51  PSD  Langley  SOL  (Autocorrelation  -  Prewhitened  and 

Smoothed) .  116 

52  PSD  Langley  CL  (Autocorrelation  -  Prewhitened  and 

Smoothed) . . .  1 1 7 

5  3  PSD  Langley  50R  (Autocorrelation  -  Prewhitened  and 

Smoothed) . 118 

54  Power  Spectral  Density  of  Altus,  75L,  Filtered 

at  500  ft . 120 

55  Powe  r  Spectral  Density  of  Altus,  Center,  Filtered 

at  500  ft .  121 

56  Power  Spectral  Density  of  Altus,  75I\,  Filtered 

at  500  ft  . 122 

57  Power  Spectral  Density  of  Bakalar,  50L,  Filtered 

at  500  ft  .  123 

58  Power  Spectral  Density  of  Bakalar,  Center,  Filtered 

at  500  ft  .  124 

59  Puwer  Spectral  Density  of  Bakalar,  50R,  Filtered 

at  500  ft  .  125 

60  Power  Spectral  Density  of  Carswell,  87 L,  Filtered 

at  500  ft . . .  126 


LIST  OF  ILLUSTRATIONS  (continued) 


IGURE 

Page 

6! 

Power  Spectral  Density  of  Carswell,  13L,  Filtered 
at  500  ft  . .  , 

127 

62 

Power  Spectral  Density  of  Carswell,  1  3R,  Filtered 
at  500  ft  . . . . 

128 

63 

Power  Spectral  Density  of  Castle,  75L,  Filtered 
at  500  ft . 

129 

64 

Power  Spectral  Density  of  Castle,  Center,  Filtered 
at  500  ft . . . 

130 

65 

Power  Spectral  Density  of  Castle,  75R,  Filtered 
at  500  ft . 

131 

66 

Power  Spectral  Density  of  Edwards,  75L,  Filtered 
at  500  ft  . 

132 

67 

Power  Spectral  Density  of  Edwards,  Center,  Filtered 
at  500  ft . 

133 

68 

Power  Spectral  Density  of  Edwards,  75R,  Filtered 
at  500  ft  . 

134 

69 

Power  Spectral  Density  of  Ellsworth,  Center,  Filtered 
at  500  ft  . 

135 

70 

Power  Spectral  Density  of  George,  50L,  Filtered 
at  500  ft  . . . . 

136 

71 

Powei  Spectral  Density  of  George,  Center,  Filtered 
at  500  ft  . . . 

137 

72 

Power  Spectral  Density  of  George,  50R,  Filtered 
at  500  ft  . .  . 

138 

73 

Power  Spectral  Density  of  Glasgow,  75L,  Filtered 
at  500  ft . . . 

139 

74 

Power  Spectral  Density  of  Glasgow,  Center,  Filtered 
at  500  ft . . . 

140 

75 

Power  Spectral  Density  of  Glasgow,  75R,  Filtered 
at  500  ft  ...............  . 

141 

76 

Power  Spectral  Density  of  Langley,  50L,  Filtered 
at  500  ft  . 

142 

77 

Power  Spectral  Density  of  Langley,  Center,  Filtered 
at  500  ft  . . . . 

143 

U 


LIST  OF  ILLUSTRATIONS  (continued) 


FIGURE 

Page 

78 

Power  Spectral  Density  of  Langley,  5QR,  Filtered 
at  500  ft . 

144 

79 

Power  Spectral  Density  of  McGuire,  75L,  Filtered 
at  500  ft  . 

145 

80 

Power  Spectral  Density  of  McGuire,  Center,  Filtered 
at  500  ft . 

146 

8! 

Power  Spectral  Density  of  McGuire,  75R,  Filtered 
at  500  ft . . . 

147 

82 

Power  Spectral  Density  of  Palmdale,  SOL,  Filtered 
at  500  ft . 

148 

83 

Power  Spectral  Density  of  Palmdale,  Center,  Filtered 
at  500  ft  . 

149 

84 

Power  Spectral  Density  of  Palmdale,  50R,  Filtered 
at  500  ft . 

150 

85 

Power  Spectral  Density  of  Shaw,  50 L,  Filtered 

at  500  ft . 

151 

86 

Power  Spectral  Density  of  Shaw,  Center,  Filtered 
at  500  ft . 

152 

87 

Power  Spectral  Density  of  Shaw,  50R,  Filtered 

at  500  ft . ...... . 

153 

88 

Power  Spectral  Density  of  Torrejon,  Center,  Filtered 
at  500  ft . . . 

154 

89 

Power  Spectral  Density  of  Torrejon,  12L,  Filtered 
at  500  ft . .  .  .  . . .  ,  . 

155 

90 

Power  Spectral  Density  of  Torrejon,  12R,  Filtered 
at  500  ft  . . . . . 

156 

91 

Power  Spectral  Density  of  Travis,  75L,  Filtered 
at  500  ft . 

157 

92 

Power  Speetral  Density  of  Travis,  Center,  Filtered 
at  500  ft  .  . . 

1 58 

93 

Power  Spectral  Density  of  Travis,  ?5L,  Filtered 
at  500  ft  .  .  . . 

159 

94 

RMS  Histories  fer  Altus  and  ftakalar  Filtered 
at  500  ft . . . . 

161 

LIST  OF  ILLUSTRATIONS  (continued) 


FIGURE  Page 

95  RMS  Histories  for  Carswell  and  Edwards  Filtered 

at  500  ft  . . .  lAd 

96  RMS  Histories  for  Ellsworth  and  Glasgow  Filtered 

at  500  ft  .  ’3 

97  RMS  Histories  for  Langley  and  McGuire  Filtered 

at  500  ft  .  ’6  4 

98  RMS  Histories  for  Shaw  and  Travis  Filtered 

at  500  ft . 165 

99  RMS  Histories  for  Torre] on  and  Castle  Filtered 

at  500  ft  . 166 

100  RMS  Histories  for  George  and  Palmdale  Filtered 

at  500  ft  .  167 

10  i  Cumulative  Distribution  of  Amplitudes  from  Altus 

Filtered  at  100  ft . 169 

I0<i  Cumulative  Distribution  of  Amplitudes  from  Bakalar 

Filtered  at  100  ft  .  170 

10  3  Cumulative  Distribution  of  Amplitudes  from  Carswell 

Filtered  at  100  ft  .  . .  171 

104  Cumulative  Distribution  of  Amplitudes  from  Edwards 

Filtered  at  100  ft . 172 

105  Cumulative  Distribution  of  Amplitudes  from  Castle 

Filtered  at  100  ft  .  . . . .  17  3 

106  Cumulative  Distribution  of  Amplitudes  from  Ellsworth 

Filtered  at  100  ft  .  17* 

107  Cumulative  Distribution  of  Amplitudes  from  George 

Filtered  at  100  ft . 170 

108  Cumulative  Distribution  of  Amplitudes  from  Glasgow 

Filtered  at  100  ft .  176 

109  Cumulative  Distribution  of  Amplitudes  from  Langley 

Filtered  at  100  ft . * . .  .  .  177 

110  Cumulative  Distribution  of  Amplitudes  from  McGuire 

Filtered  at  100  ft . 178 

111  Cumulative  Distribution  of  Amplitudes  from  Palmdale 

Filtered  at  100  ft  .  179 


LIST  OF  ILLUSTRATIONS  (concluded) 


FIGURE  Page 

112  Cumulative  Distribution  of  Amplitudes  from  Shaw 

Filtered  at  100  ft  .................  180 

]  13  Cumulative  Distribution  of  Amplitudes  from  Torrejon 

Filtered  at  100  ft  . .  181 

114  Cumulative  Distribution  of  Amplitudes  from  Travis 

Filtered  at  100  ft  . .  182 

115  RMS  Histories  for  Altus  and  Bakalar  Filtered 

at  100  ft . 184 

116  RMS  Histories  for  Carswell  and  Edwards  Filtered 

at  100  ft  . 185 

117  RMS  Histories  for  Ellsworth  and  Glasgow  Filtered 

at  100  ft  . . 186 

118  RMS  Histories  for  Langley  and  McGuire  Filtered 

at  100  ft  . . 187 

119  RMS  Histories  for  Shaw  and  Travis  Filtered 

at  100  ft .  188 

120  RMS  Histories  for  Torrejon  and  Castles  Filtered 

at  100  ft . . . .  .  189 

121  RMS  Histories  for  George  and  Palmdale  Filtered 

at  100  ft . 190 


LIST  OF  TABLES 


TABLE 

Page 

I 

Identification  of  Data  Analyzed  . 

7 

II 

Filter  Constants  . . . 

10 

III 

Confidence  Limit  Factors  for  Spectral  Estimates  .  .  . 

12 

IV 

RMS  Values  for  Data  Filtered  at  500  ft  . 

25 

V 

Number  of  Profiles  with  Relative  Maxima  of  PSD 
at  Indicated  Ranges  of  Wavelength  .......... 

31 

VI 

Mean  Square  Values  for  Segments  of  Profiles 

Filtered  at  500  ft  . . . 

39 

VII 

RMS  Values  for  Data  Filtered  at  100  ft . .  . 

44 

VIII 

Cumulative  Distributions  of  Amplitudes  for  Half 
Wavelength  Intervals  .  . . . 

53 

XI 

Distribution  of  Long  Wavelength  Peaks  Input 
to  Simulated  Runway  . . 

57 

X 

Distribution  of  Long  Wavelength  Peaks  Output 

from  Simulated  Runway . 

57 

XI 

Simulated  Runway  Elevations  in  Inches  at 

Two-Foot  Intervals  . . 

63 

XII 

Filter  Guidelines  . . . 

84 

xiii 


LIST  OF  SYMBOLS 


Symbol 


D.  F. 


FREQ 


P{A} 


Description 


Distance  down  runway  (ft) 

Degrees  of  freedom 

Cutoff  frequency  in  low-pass  filter 

Nyquist  frequency 

Sampling  rate  (points/ft) 

Amplitude  as  function  of  increment 
down  runway  (in) 

Fourier  transform  of  x(d) 

Fast  Fourier  transform 

F  requency 

Filter  weight 

Fourier  transform  of  filter  weights 

Number  of  data  values  in  transform 
when  computing  PSD  or  the  number 
of  filter  weight  values 

Probability  of  A 

Power  spectral  density 

Root  mean  square 

Radian 

Original  runway  elevations  as 
function  of  distance 

Low- pass  filtered  data 

High- pass  filtered  data 

Increment  between  samples  (ft) 

Roll-off  frequency  in  low-pa3s 
filter 

Estimate  of  filter  error 
Wavelength  (ft) 

Cutoff  wavelength  in  high-pass 
filter 


i 

't  k 

:'i  % 

i  H 

1  1 

i  I 


|  -§j 

,J  ;5j 

LIST  OF  SYMBOLS  (continued)  *  1 


I 


Symbol 

Description 

O’ 

Standard  deviation 

«(0) 

Power  spectral  density  as  function 
of  reduced  frequency  (in^/rad  ft) 

$*(£2) 

Normalized  power  spectral  density 
as  function  of  reduced  frequency 

4>(“) 

Power  spectral  density  as  function 
of  frequency  (cycles/ft) 

.Lt 

i  1 
•| 

CO 

Frequency  (cycles /ft) 

J2 

Reduced  frequency  (rad/ft) 

?  1 

SECTION  1 
INTRODUCTION 


Runway  roughness  has  long  been  recognized  as  a  source  of  dynamic 
excitation  for  an  aircraft.  Changes  in  vertical  runway  profiles  have  been 
identified  not  only  as  major  sources  of  fatigue  damage  for  some  aircraft 
types,  but  also  as  a  potential  cause  of  loss  of  control  due  to  magnification  of 
forces  as  they  are  transmitted  through  the  aircraft  structure  to  the  cockpit. 
The  effect  of  a  runway-induced  excitation  obviously  depends  on  the  runway 
profile,  the  properties  of  the  landing  gear,  and  the  dynamic  characteristics 
of  the  aircraft. 

Typically,  at  some  stage  in  the  development  of  an  aircraft,  a  mathe¬ 
matical  model  is  available  to  a  designer  or  structural  analyst  which,  when 
coupled  with  a  given  runway  profile,  is  capable  of  providing  a  description 
of  a  response  of  interest.  Thus,  given  the  model  and  the  runway  profile, 
responses  of  interest  can  be  predicted.  Further,  it  can  be  assumed  that 
the  model  can  be  transformed  so  that  calculations  can  be  performed  in 
either  the  frequency  or  time  domains.  If  these  calculations  are  being 
performed  in  the  design  stages,  changes  in  the  landing  gear  are  reflected 
by  changes  in  the  parameters  or  mathematical  terms  of  the  model.  The 
effect  of  the  change  for  the  given  runway  can  then  be  determined.  However, 
the  choice  of  the  runway  profile  or  profiles  and  the  form  of  its  mathematical 
description  that  should  be  used  in  the  exercise  of  the  model  are  quite  difficult. 
No  aircraft  type  will  operate  on  only  one  runway  and  the  roughness 
characteristics  between  runways  exhibit  significant  variation  that  is 
important  to  aircraft  response. 

There  are  three  dominant  methods  for  analyzing  the  effectiveness  of  a 
landing  gear  system  to  runway  excitation  and  these  can  be  described  as  the 
discrete,  statistical,  and  time  history  approaches.  In  the  simplest  of  these 
analyses,  the  discrete  approach,  the  runway  is  characterized  as  consisting 
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of  a  series  of  "  1  -  cos"  bumps  and  depressions  where  the  amplitudes  and 
wavelengths  of  the  disturbances  are  specified  from  an  observed  joint  distri¬ 
bution  of  bump  heights,  depression  depths  and  their  corresponding  wavelengths 
as  determined  from  existing  runway  profiles.  The  aircraft  response  of 
interest  is  calculated  for  all  relevant  combinations  of  amplitude,  wavelength 
and  aircraft  ground  velocity.  Since,  in  this  analysis,  each  combination  of 
amplitude  and  wavelength  is  treated  separately,  the  possibility  of  existing 
dynamic  conditions  at  the  time  of  entry  into  the  excitation  is  ignored.  Further, 
the  "1  -  cos"  shape  of  the  bumps  and  depressions  is  not  considered  to  be 
representative  of  actual  runways  and  there  exists  considerable  subjectivity 
in  the  determination  of  the  joint  distributions  of  amplitudes  and  wavelengths,. 

For  these  reasons,  the  results  of  a  discrete  analysis  are  considered  only  as 
a  preliminary  indication  of  the  adequacy  of  a  landing  gear  system. 

The  statistical  approach  is  centered  on  the  concept  that  a  runway  can  be 
considered  as  a  stationary,  gaussian  excitation.  In  this  mathematical 
framework,  the  characterization  of  the  runways  is  performed  in  the  frequency 
domain  by  modeling  the  power  spectral  densities  (PSD)  of  available  runway 
elevation  profiles.  Knowledge  of  the  transfer  function  of  the  random  dynamic 
excitation  to  a  response  of  interest,  coupled  with  known  theory  of  time  history 
analysis  is  then  used  to  predict  level  crossing  or  peak  distributions  of  the 
response  parameter.  This  approach  provides  a  more  realistic  prediction  of 
response  than  does  the  discrete  approach  since  the  PSD  and  transfer  function 
account  for  the  continuously  changing  excitation.  However,  it  has  been  found 
that  the  predicted  response  from  the  statistical  approach  is  generally  more 
severe  than  the  actual  when  a  given  runway  is  analyzed.  This  lack  of  agreement 
is  due  to  a  violation  of  the  assumptions  of  stationarity  or  gaussianness  or  to 
the  nonlinearity  of  the  system  which  is  also  a  required  assumption  in  the 
statistical  approach. 
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The  best  method  of  predicting  response  to  runway  excitation  is  the 
time  history  approach.  In  this  method  a  mathematical  model  relating  a 
response  of  interest  to  the  changes  in  elevation  profiles  is  derived  and  the 
actual  runway  profiles  are  used  as  the  input  to  the  model.  The  time  history 
of  the  response  is  then  calculated  and  the  accuracy  of  this  prediction  is 
dependent  strictly  on  the  mathematical  model  of  the  aircraft.  The  difficulty 
with  this  approach  is  that  an  actual  runway  profile  is  required  as  an  input. 

If  the  response  of  a  particular  aircraft  to  a  particular  runway  is  desired, 
there  is,  of  course,  no  problem.  The  elevation  profile  of  the  runway 
provides  the  input  of  interest.  However.-  in  the  design  stages  of  an  aircraft, 
the  choice  of  the  runway  profile  that  should  be  used  is  not  clear  since 
aircraft  types  respond  differently  to  a  specific  runway.  Currently,  this 
problem  is  approached  by  specifying  that  the  runway  profile  used  in  analysis 
must  have  a  power  spectral  density  which  exceeds  that  given  in  specifications. 
Most  runway  profiles  can  be  used  in  this  procedure  since  by  appropriately 
scaling  the  elevation  measurements  of  a  profile,  the  PSD  can  be  raised  or 
lowered  to  the  specified  limits.  It  can  easily  be  seen,  however,  that  this 
procedure  can  lead  to  an  overdesign.  It  may  well  happen  that  to  raise  the 
power  spectral  density  to  the  acceptable  level  at  a  noncritical  wavelength, 
the  power  at  a  critical  wavelength  may  be  greater  than  can  realistically  be 
expected. 

In  specifying  runway  roughness  design  criteria,  the  above  discussion 
leads  to  a  dilemna.  While  runways  can  be  characterized  in  terms  of  joint 
distributions  of  bump  heights,  depression  depths  and  corresponding  wave¬ 
lengths  and  also  in  terms  of  power  spectral  densities,  the  methods  of 
calculating  response  using  these  characterizations  are  unsatisfactory.  On 
the  other  hand,  the  time  history  approach  which  provides  the  best  method 
of  predicting  response  requires  a  runway  profile  as  input  and  the  choice  of 
the  appropriate  profile  is  not  clear.  This  study  represents  an  attempt  to 
resolve  the  dilemna  by  achieving  the  following  two  objectives: 
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1)  To  determine  and  evaluate  quantitative  characteristics  of  runway 
profile  elevations  which  are  pertinent  to  aircraft  ground  loads 
analyses,  and 

2)  To  develop  a  method  for  synthesizing  artificial  runways  which 
reflect  the  quantified  characterization  of  the  real  runways  which 
were  used  for  analysis. 

To  achieve  these  objectives,  4-0  lines  of  profile  elevations  from  14 
runways  were  analyzed.  To  make  the  profiles  comparable  and  to  remove 
grade  effects  of  no  dynamic  interest,  initially  all  lines  were  digitally,  high- 
pass  filtered  at  a  cutoff  wavelength  of  500  ft.  Power  spectral  densities  and 
amplitude  probability  densities  were  obtained  for  each  line.  In  order  to 
investigate  stationarity,  power  spectral  densities  were  calculated  for  each 
third  of  each  runway.  Since,  for  some  of  the  runways,  the  power  spectral 
densities  of  the  segments  were  significantly  different,  further  study  of  the 
lack  of  stationarity  was  investigated  by  calculating  the  RMS  value  for  each 
512  ft  section  each  runway.  These  RMS  values  also  displayed  significant 
variation  and  it  was  determined  that,  generally,  the  significantly  larger 
values  were  the  result  of  long  wavelength  undulations  in  sections  of  the 
runways.  It  was  concluded,  therefore,  that  stationarity  could  not  be  assumed 
for  the  profiles  filtered  at  500  ft  since  they  were  not  self- stationary. 

Further,  the  investigation  of  the  amplitude  probability  densities  could  not 
generally  be  modeled  by  a  gaussian  distribution  primarily  because  of  an 
excess  of  values  in  the  tails  of  the  distribution.  This  excess  of  probability 
in  the  tails  would  also  be  a  result  of  the  long  wavelength  undulations. 

The  profiles  were  then  high-pass  filtered  at  a  cutoff  wavelength  of 
100  ft  and  it  was  concluded  that  for  wavelengths  less  than  100  ft  the 
assumptions  of  stationarity  and  gaussianness  were  far  more  reasonable. 
Therefore,  a  compound  characterization  was  derived  in  which  the  wave¬ 
lengths  less  than  100  ft  were  modeled  by  the  PSD's,  assuming  stationarity 
and  gaussianness,  and  the  wavelengths  between  100  ft  and  500  ft  were 
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modeled  in  a  discrete  fashion  by  determining  the  joint  distributions  of  primary 
peaks,  associated  wavelength,  and  location  on  the  runway  from  the  profiles 
filtered  at  a  cutoff  wavelength  of  500  ft.  This  compound  characterization  was 
then  used  to  synthesize  a  simulated  runway. 

The  details  and  results  of  these  operations  arc  presented  in  the 
following.  First,  the  computational  and  data  processing  techniques  are 
presented  in  Section  2.  The  results  of  the  analysis  of  the  data  filtered  at 
500  ft  and  at  100  ft  are  presented  in  Section  3,  And  the  method  of  synthesizing 
the  simulated  runway  and  a  sample  result  are  presented  in  Section  4.  A 
brief  summary  and  conclusions  are  contained  in  Section  5. 
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SECTION  2 
DATA  PROCESSING 


The  following  is  a  description  of  the  procedures  used  in  the  processing 
and  analysis  of  the  runway  roughness  data.  Appendix  I  gives  a  more  detailed 
description  of  the  digital  filtering  techniques.  All  processing  of  the  data  was 
carried  out  on  the  CDC  6600  computer  located  at  Wright-Patterson  AFB. 

2.  1  DESCRIPTION  OF  RAW  DATA 

The  raw  data  as  received  was  in  the  form  of  runway  elevations  at 
intervals  of  either  six  inches  or  two  feet.  The  data  for  Torrejon  is  the 
only  data  at  two-foot  intervals  as  measured  by  rod  and  level  and  punched 
on  cards.  All  of  the  other  data  was  measured  and  recorded  at  six-inch 
intervals  by  an  automatic  runway  profile  measuring  instrument,  profilometer, 
on  magnetic  tape  and  later  punched  on  cards.  The  rod  and  level  elevations 
were  in  feet,  which  were  subsequently  converted  to  inches,  and  the 
profilometer  elevations  in  inches.  All  elevations  are  relative  to  some 
arbitrary  starting  point. 

The  description  of  the  lines  of  survey  and  the  total  number  of  data 
points  for  each  line  is  given  in  Table  I. 

2.  2  FILTERING  METHODS 

The  data  from  all  of  the  runways  contained  a  linear  trend  which 
biases  the  power  spectral  density  considerably.  This  linear  trend  must  be 
removed  and  digital  filtering  was  considered  the  most  reliable  method  to 
accomplish  this.  The  digital  filter  that  was  designed,  consisted  of  a  two 
stage  Martin-Graham  low-pass  filter,  the  results  of  which  were  subtracted 
from  the  original  data  so  as  to  provide  high-pass  filtering  of  the  data,  A 
two  stage  filter  was  needed  to  provide  a  relatively  short  roll-off.  The 
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TABLE I 

IDENTIFICATION  OF  DATA  ANALYZED 


Altus 


Bakalar 


Carswell 


Edwards 


Ellsworth 

Glasgow 


Langley 


McGuire 


Shaw 


Travis 


Torrejon 


Castle 


George 


Palmdale 


Runway 


•Alter  interpolation  to  6  inch  sete^'^le 


No.  Points 

26,  930 

26,  910 

25,  190 

10,  000 
10, 000 
9,  600 

24,  060 
24,  100 
24,  060 

28,  800 
28,  800 
28,  800 

27,  120 

26.  960 
27,010 
15,030 

20,  040 
20,050 
19,010 

14,  960 
14.  950 
15,070 

20. 070 
18, 880 
18. 000 

22.  080 
22, 040 

22.  040 

26. 790* 
26, 710* 
26. 790* 

23,  650 
23,610 
23. 660 

18.  260 
16,  270 
18, 000 

24,030 

24,030 

24,03) 


.■aiwwn 


filtering  was  performed  with  the  following  equation  which  is  developed  in 
Appendix  I. 


*2  "l 

^  ^  ^1  X{d+ni  Adi+n2AdA 

n2=-N2  *  nj-'-Nj  1  * 


The  filter  constants  used  are  shown  in  Table  II.  The  response  of  the  high- 
pass  filter  is  shown  in  Figure  1.  The  500  ft  wavelength  cutoff  was  chosen 
because  it  was  believed  that  greater  wavelengths  were  ot  no  interest.  The 
roll-off  was  accomplished  in  as  short  a  frequency  span  as  practical.  The 
data  for  each  line  of  survey  was  filtered  and  plotted.  Upon  scanning  those 
plots,  numerous  single  point  spikes  were  noticed.  These  were  edited  from 
the  filtered  data  on  the  premise  that  they  could  not  be  real  data  if  the 
amplitude  was  greater  than  0.4  inches. 


After  some  further  analysis  it  was  decided  that  the  filter  cutoff 
frequency  would  be  increased  to  correspond  to  a  wavolength  of  100  ft. 
Filtering  of  all  lines  of  survey  was  then  accomplished  using  Equation  (1)  and 
constants  in  Table  II  which  provide  a  high-pass  filter  with  the  response 
shown  in  Figure  2. 


2.  3  PSD  COMPUTATION 

The  Power  Spectral  Density  (PSD)  of  each  line  of  survey  was  computed 
using  the  Fast  Fourier  Transform  (FFT)  on  the  filtered  data.  In  order  to 
get  more  degrees  of  freedom  of  the  PSD  estimates,  the  PSD  values  were 
computed  for  each  512  ft  segment  of  a  line  of  survey  and  these  estimates 
averaged  to  obtain  the  PSD  for  the  entire  line.  In  order  to  make  the  most 
efficient  use  of  the  FFT  to  compute  the  PSD  coefficients,  it  is  necessary 
that  the  data  array  siae  be  an  integer  power  of  2,  The  array  sine  chosen 
was  1024,  which  explains  the  previously  mentioned  512  ft*  The  sine  of 
this  array  also  sets  the  frequency  resolution  of  the  PSD  output  data; 
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Response  of  500  ft  Cutoff  Filter  Figure  2  Response  of  100  ft  Cutoff  Filter 


TABLE  II 

FILTER  CONSTANTS 


. 

500  ft  Cutoff 

_ _ _ __ _ _ 

100  ft  Cutoff 

Preliminary 

Final 

Preliminary 

Final 

f. 

2 

0.0667 

2 

0.2 

Ad 

0.5 

15 

0.5 

5 

fc 

0.0015 

0.0015 

0.007 

0.  007 

Af 

0.030 

0.0005 

0.033 

0. 003 

N 

75 

140 

80 

150 

A FREQ  = 


sample  rate  _  2 

array  sisse  ~  1024 


0.001953  cycles/foot  or 
0.01227  rad/ foot 


The  PSD  coefficients  were  computed  using  the  following  equation: 


«(^j 


a 


2w(N-l) 


.  (2) 


The  2  in  the  numerator  of  Equation  >;2)  is  required  to  distribute  the  total 
power  over  the  positive  frequencies  only  and  the  2v  in  the  denominator  is  to 
convert  the  value  to  reduced  frequency.  The  units  of  the  PSD  values  are  in 
in *•/( rad/ft).  The  number  of  degrees  of  freedom  for  each  of  the  PSD's  is 
equal  t->  the  number  of  512  ft  sections  in  the  data  times  two.  This  number 
can  vat  sv  two  on  a  given  pint  as  a  function  of  frequency  since  a  PSD  was 
computed  for  that  portion  of  a  line  leftover  after  some  integral  multiple  of 
SI 2  ft  sections  were  used.  For  this  last  piece  of  data  on  each  Une*  only 
{he  P$D‘s  for  frequencies  having  at  least  one  full  cycle  were  o^ed  in  the 
averaging  process.  Therefore,  a  given  line  can  have  2N  degrees  of  freedom 
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up  to  some  frequency  and  2(N  +  1)  degrees  of  freedom  for  higher  frequency, 
where  N  is  the  number  of  integral  multiples  of  512  ft. 

By  breaking  the  line  into  512  ft  sections,  the  PSD's  of  various  sections 
of  the  runway  could  be  computed.  Each  of  the  lines  was  broken  up  into 
one-thirds  and  the  PSD's  computed  for  each  third  to  determine  if  the  ends 
of  the  runways  were  rougher  than  the  centers.  In  order  to  have  an  equal 
number  of  degrees  of  freedom  for  each  of  the  thirds,  it  was  sometimes 
necessary  to  include  adjacent  512  ft  sections  in  two  of  the  thirds. 

To  clarify  any  confusion  that  might  arise  on  the  use  of  the  FFT  to 
compute  PSD's,  a  comparison  was  made  between  the  PSD's  computed  by 
talcing  the  transform  of  the  autocovariance  function  with  the  PSD's  obtained 
by  use  of  the  FFT.  Since  there  was  a  program  available  to  compute  the 
PSD  using  the  autocovariance  function  which  also  contained  prewhitening 
and  smoothing,  the  FFT  method  was  reprogramed  to  include  prewhitening 
and  smoothing  for  comparison  purposes.  The  results  of  this  comparison 
are  shown  in  Appendix  II,  The  number  of  degrees  of  freedom  associated 
with  each  PSD  is  shown  on  the  plot.  Confidence  limits  for  the  PSD  plots 
are  a  function  of  the  number  of  degrees  of  freedom  on  the  estimates  and 
can  be  determined  from  Table  III,  Using  this  table,  it  can  be  seen  that  for 
50  degrees  of  freedom  there  is  80%  confidence  that  the  true  value  falls 
within  0.  75  times  the  estimate  and  1.  26  times  the  estimate. 

2.4  LONG  WAVELENGTH  DESCRIPTION 

The  term  "long  wavelength"  is  used  here  to  describe  all  primary  peaks 
and  valleys  greater  than  100  ft  wavelength.  These  peaks  and  valleys  were 
extracted  from  the  data  after  it  had  been  filtered  at  the  500  ft  cutoff.  The 
method  used  selects  the  maximum  excursion  between  crossings  of  the  zero 
mean  line.  Only  peaks  with  an  amplitude  greater  than  0.  25  inches  were 
retained.  These  peaks  were  assigned  a  wavelength  equal  to  twice  the 
distance  between  mean  crossings.  Peaks  were  then  tabulated  by  half 


11 


TABLE  III 

CONFIDENCE  LIMIT  FACTORS  FOR  SPECTRAL  ESTIMATES 


Degrees  of 

F  reedom 

Exceeded  by  90% 
of  all  values 

Exceeded  by  50% 
of  all  values 

Exceeded  by  10% 
of  all  values 

10 

0.  49 

0.  93 

1.60 

20 

0.  62 

0.96 

1.42 

30 

0.69 

0.  98 

1.  34 

40 

0.  73 

0.  98 

1.  30 

50 

0.  75 

0.99 

1.26 

100 

0.  82 

0.  99 

1.  18 

wavelength  and  amplitude  into  two  tables,  one  containing  positive  peaks  and 
the  other  containing  negative  peaks.  In  the  process  of  selecting  the  peaks, 
all  wavelengths  were  retained. 

It  is  thought  that  using  primary  peaks  to  define  the  long  wavelengths 
in  conjunction  with  the  PSD  representation  of  the  shorter  wavelengths  is 
the  most  logical  method  of  representing  the  runway  profiles.  It  is  believed 
that  most  of  the  data  greater  than  100  ft  wavelength  is  in  the  form  of  primary 
peaks.  An  earlier  attempt  at  selecting  cyclic  information  from  the  data 
produced  relatively  large  wavelengths  that  were  known  to  have  been  removed 
from  the  data. 

2.  5  OTHER  DESCRIPTORS 

Amplitude  probability  densities  were  computed  for  all  lines  of 
survey.  These  were  computed  in  the  same  program  that  computed  PSD. 

Root  mean  square  (RMS)  values  were  computed  for  each  line  of 
survey  and  for  each  512  ft  section  of  each  line.  The  RMS  values  for  each 
entire  line  of  survey  were  computed  in  two  ways.  First,  the  RMS  value 


12 


was  computed  on  a  point  for  point  basis  from  the  filtered  data  and  was  also 
computed  as  the  square  root  of  the  area  under  the  PSD  curve.  These  agree 
to  two  significant  figures  in  all  cases.  The  RMS  values  for  the  512  ft 
sections  and  for  the  profiles  broken  into  thirds  were  computed  from  the 
square  root  of  the  area  under  the  PSD  curves  only. 
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SECTION  3 
DATA  ANALYSIS 


The  basic  rationale  in  the  characterization  of  the  available  runway 
elevation  profiles  is  that  these  data  are  a  representative  sample  of  all 
concrete  runways  which  will  be  routinely  traversed  by  Air  Force  aircraft 
during  their  lifetime.  Any  one  runway,  of  course,  provides  a  reasonable 
fixed  dynamic  excitation  to  an  aircraft  in  that  its  elevation  profile  would  be 
expected  to  change  only  slowly  with  time.  But  since  at  the  design  stages  of 
an  aircraft,  it  is  not  known  which  of  the  multitude  of  possible  ground 
surfaces  the  aircraft  will  traverse,  a  statistical  characterization  of  runway 
profiles  becomes  necessary.  Thus,  although  all  of  the  available  elevation 
profiles  have  been  resuriaced  since  the  data  were  collected  and  although 
multiple  lines  wore  obtained  from  the  same  runways,  it  is  assumed  in  the 
analysis  that  each  line  of  data  is  representative  of  the  population  of  possible 
inputs  to  an  aircraft  during  takeoff  and  landing. 

The  original  elevation  measurements  contain  the  long  term  grades 
exhibited  by  the  runways  (c,  f.  Reference  1).  Assuming  that  any  grade 
which  has  a  wavelength  greater  than  500  ft  would  have  no  dynamic  input  to 
a  modern  aircraft,  the  long  term  grades  for  all  lines  of  data  were  removed 
by  a  high-pass  digital  filter  which  essentially  passed  all  wavelengths  (\) 
less  than  500  ft  and  removed  all  wavelengths  greater  than  666  ft.  Between 
500  ft  and  666  ft,  the  amplitudes  were  attenuated  as  described  in 
Paragraph  2.  2.  The  filtered  profiles  provided  the  basic  data  set  which 
were  then  analysed  with  particular  emphasis  on  the  statistical  assumptions 
required  for  the  power  spectral  density  characterization.  The  results  of 
this  analysis  are  presented  in  paragraph  3. 1,  Since  it  was  found  that,  in 
general,  the  profiles  filtered  at  a  500  ft  are  not  stationary  or 

gaussian,  the  profiles  were  again  filtered  with  *maJS  3  100  ft.  The  analysis 
of  the  profiles  with  the  lower  cutoff  wavelength  is  presented  in  Paragraph  3.  2 


14 


and  the  discrete  characterization  of  the  wavelengths  between  100  and  500  ft 
is  presented  in  Paragraph  3.  3. 

3,  1  MAXIMUM  WAVELENGTH  OF  500  FEET 

In  determining  the  significant  characteristics  of  runway  roughness  it 
is  apparent  that  the  most  desirable  characterization  is  that  which  requires 
the  fewest  parameters  while  retaining  sufficient  information  for  accurate 
ground  loads  analysis.  Since  from  the  response  viewpoint  the  important 
roughness  characteristics  are  the  magnitude  and  corresponding  wavelength 
of  the  surface  undulations,  the  power  spectral  density  characterization  of 
the  runways  appears  to  provide  a  workable  model.  The  rationale  of  this 
model  is  that  each  individual  runway  is  a  representation  of  a  stationary, 
gaussian  input  and  the  family  of  power  spectral  densities  can  be  modeled  as 
a  function  of  the  RMS  values  of  the  elevation  profiles  and,  perhaps,  two  cr 
three  othor  parameters.  Past  studies  have  indicated  a  similarity  in  shape 
of  the  power  spectral  densities  from  runway  to  runway,  making  such  a 
model  at  least  plausible.  However,  when  the  statistical  response  predicted 
by  the  power  spectral  density  of  a  runway  is  compared  to  the  response 
predicted  by  a  time  history  analysis,  the  power  spectral  density  generally 
overpredicta  the  magnitude  of  loads.  Therefore,  in  the  initial  stages  of 
this  study  considerable  emphasis  was  placed  on  the  power  spectral  density 
characterization  of  runway  roughness  with  particular  attention  on  the 
assumptions  required  by  this  approach. 

Removing  the  grades  present  in  the  raw  data  by  digitally  filtering  at 
the  longest  wavelength  of  interest  {500  ft)  results  in  an  elevation  profile 
of  essentially  aero  mean.  This  permits  comparison  of  amplitude 
probability  densities  across  runways  and  allows  a  preliminary  investigation 
of  statioaarity  by  determining  if  the  RMS  values  of  sections  of  one  runway 
can  be  considered  as  constant.  The  following  paragraphs  are  devoted  to 
an  analysis  of  the  amplitude  probability  density  function,  the  power  spectral 


densities  of  the  profiles,  and  the  study  of  stationarity  for  the  profiles 
filtered  at  500  ft. 

3.  1.  1  Amplitude  Probability  Densities 

An  important  element  in  the  characterization  of  any  time 
series  is  the  distribution  of  the  amplitude  values  about  their  mean  value.  In 
considering  a  runway  profile  to  be  a  representative  sample  of  random  time 
series,  the  assumption  of  ergodicity  is  required.  This  assumption  permits 
the  estimation  of  statistical  properties  from  a  single  time  history  as 
opposed  to  obtaining  the  estimates  from  an  ensemble  of  time  histories. 

The  amplitude  probability  density  can  be  obtained  quite  simply  for  any 
profile  by  counting  the  number  of  data  points  in  defined  ranges  of  deviation 
from  the  mean.  The  distribution  thus  obtained  describes  the  amplitudes 
without  regard  to  wavelength  and  can  also  be  used  to  obtain  the  RMS  value 
and  other  moments  of  the  time  history.  Further,  most  of  the  available 
mathematics  for  the  application  of  time  theory  methods  to  practical 
problems  depends  on  the  assumption  that  the  amplitude  probability  density 
is  gaussian. 


In  considering  the  amplitude  probability  densities  of  runway 
profiles  it  is  easily  seen  that  the  definition  of  the  mean  value  is  critically 
important.  For  example,  if  grades  present  in  the  original  profiles  were 
not  removed,  the  amplitude  probability  densities  for  some  of  the  profiles 
would  appear  as  a  uniform  distribution.  Runway  35  at  Altus  AFB  has  an 
approximately  linear  grade  of  350  inches  over  the  length  of  the  runway 
(c.f.  Reference  1).  The  relatively  smaller  amplitudes  which  are  important 
to  aircraft  would  be  completely  lost  if  this  grade  (or  nonrandom  component) 

a. 

were  not  removed.  Thus,  if  x^  represents  the  original  profile,  x'^ 
represents  the  nonrandom  or  nonrelevant  portion  of  the  profile,  the  reLevant 
portion  is  defined  by 


=  *d  - 


(3) 


i 
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That  is,  x'^  is  the  mean  value  about  which  the  amplitude  deviations  are 
considered  as  important  to  roughness  characterizations.  As  noted  previously, 
this  paragraph  is  devoted  to  the  analysis  of  amplitudes  about  a  mean  defined 
as  a  result  of  a  low-pass  filtering  operation  in  which  all  wavelengths  less 
than  500  ft  were  removed. 

Of  prime  concern  in  the  investigation  of  these  amplitudes  is 
the  applicability  of  one  family  of  distributions  (e.  g. ,  gaussian).  That  is, 
can  one  family  be  used  to  model  the  amplitude  probability  densities  of  all 
the  profiles  with  the  differences  in  the  profiles  being  reflected  by  the  values 
of  the  parameters  of  the  family.  Particular  attention  was  given  to  the 
gaussian  distribution  because  of  the  wealth  of  available  theory  for  this 
distribution  and  because  of  the  simplicity  of  synthesizing  a  time  history 
with  any  given  power  spectral  density  when  the  amplitude  probability 
densities  are  gaussian. 


The- amplitude  probability  densities  for  all  lines  of  data  were 
converted  to  cumulative  distributions  and  plotted  on  normal  probability 
paper.  On  this  particular  type  of  paper  an  exact  gaussian  distribution  is 
represented  by  a  straight  line  with  the  slope  of  the  line  being  a  function  of 
the  standard  deviation  and  the  abscissa  of  the  50th  percentile  being  the 
mean.  In  general,  the  mid  ranges  of  the  distributions  (between  the  5th  and 
95th  percentiles)  were  linear  but  deviations  in  tails  were  usually  present. 

To  provide  an  indication  of  the  types  of  distributions  observed.  Figures  3, 

4  and  5  present  the  cumulative  distributions  on  normal  probability  paper 
for  six  of  the  40  lines  of. data.  The  solid  lines  on  those  plots  are  those 
derived  from  the  calculated  means  and  standard  deviations  of  the  amplitudes 
as  opposed  to  the  best  straight  edge  fit  of  data  points. 


---  •- 
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Percent  of  data  points  exceeding  indicated  amplitude 


Percent  of  data  points  exceeding  Indicated  amplitude 


The  distributions  of  Figure  3  indicate  that  the  gaussian 
distribution  provides  an  excellent  model  for  the  amplitude  probability 
densities  of  these  two  lines  of  data.  It  should  be  noted  there  is  no 
statistical  test  of  the  goodness  of  fit  of  these  distributions  since  the 
amplitudes  are  not  statistically  independent.  However,  in  view  of  the 
large  number  of  data  points  involved,  close  agreement  would  be  required 
for  an  acceptable  fit.  This  degree  of  agreement  between  the  sample 
distributions  and  the  gaussian  model,  however,  was  uncommon  among  the 
40  lines  of  data.  Figure  4  displays  the  cumulative  distributions  of  two 
lines  for  which  the  gaussian  model  provides  an  inadequate  fit  of  the 
amplitude  probability  densities.  Note  in  this  figure  that  the  gaussian 
distribution  displays  a  marked  disagreement  at  the  mid  values  as  well  as 
in  the  tails.  It  may  bo  of  interest  to  note  that  one  of  the  lines  from  Travis 
closely  follows  a  gaussian  distribution  while  another  Line  from  the  same 
runway  is  not  gaussian.  The  linos  displayed  in  Figure  5  are  considered 
representative  of  the  majority  of  the  amplitude  distributions.  They  display 
reasonable  agreement  in  the  tnid  ranges  and  disagreement  in  the  tails.  The 
distribution  for  the  runway  from  Bakular  is  slightly  skewed  to  the  right 
while  the  distribution  for  the  runway  from  Torrejon  is  slightly  skewed  to 
tile  left. 

These  figures  display  clearly  that  no  single  family  of  distribu¬ 
tions  would  provide  a  good  model  for  all  of  the  available  line*  of  data.  They 
are  both  positively  and  negatively  skewed  indicating  a  change  in  the  shape 
parameter  of  a  distribution  and  although  the  gaussian  model  provides  a 
goad  fit  to  a  few  of  the  lines,  it  does  not  provide  a  good  fit  for  others. 
Therefore,  at  least  two  models  would  he  required  to  describe  all  of  the 
lines  satisfactorily  and  the  complexity  introduced  by  multiple  models  is 
sufficiently  great  that  fitting  each  line  is  unwarranted  at  present.  If  these 
amplitude  probability  densities  must  be  modeled,  >:he  gaussian  distribution 
would  provide  as  good  a  fit  as  any  other  distributional  model  and,  considering 


the  available  theory  for  time  history  analysis  for  the  gaussian  model,  would 
be  recommended. 

One  other  point  regarding  the  amplitude  probability  densities 
should  be  noted  at  this  time.  The  elevation  profiles  of  the  lines  which  are 
not  gaussiau  each  contain  several  cycles  of  high  amplitude,  long  wave¬ 
length  undulations  that  are  not  typical  of  the  other  sections  of  the  runways. 
Conversely,  the  elevation  profiles  of  the  lines  which  have  gaussian  amplitude 
probability  densities  do  not  contain  such  cycles.  This  observation  is  taken 
as  an  indication  that  high  amplitude,  long  wavelength  cycles  which  are  not 
consistently  presont  in  a  runway  profile  tends  to  cause  nongaussiatmess  in 
the  amplitude  probability  density.  Tnis  contention  is  further  supported  by 
the  realization  that  if,  say,  20  percent  of  a  runway  profile  consisted  of 
high  amplitude,  long  wavelength  cycles  while  the  remaining  80  percent 
consisted  of  a  homogeneous  mixture  of  lower  amplitudes  and  wavelengths, 
the  amplitude  probability  density  would  be  distorted  by  the  high  values. 
Further,  the  RMS  value  for  the  entire  runway  would  be  representative  of 
neither  type  of  amplitudes.  Therefore,  the  somewhat  typical  pattern 
observed  in  the  amplitude  distributions  can  be  explained  by  the  presence 
of  high  amplitude,  long  wavelength  undulations  of  the  profiles.  In 
Paragraph  5.  2,  1,  the  amplitude  probability  densities  ol  profiles  filtered 
at  a  cutoff  wavelength  of  ‘>;Q  ft  are  discussed  and  it  is  shown  that  the 
gaussian  distribution  provides  an  acceptable  model  for  these  densities. 


Since  the  RMS  values  for  the  elevation  profiles  provide  a 
gross  measure  of  the  runway  roughness,  these  values  are  presented  in 
Table  IV  for  each  of  the  40  lines  of  data.  The  histogram  of  the  RMS 
values  is  presented  in  Figure  6  to  display  the  distributional  properties  of 
this  parameter.  For  this  distribution  the  moment  valuos  are  given  by 

mean  =  0.  3f>4 

standard  deviation  =  0.  176 

skewness  =  2.  306 

kurtosis  =  3.  991 

A  Kolmogroff-Smirnov  goodness  of  fit  test  rejects  the  hypothesis  that  the 
distribution  is  qaur.sian  at  a  99  percent  level  of  confidence.  Therefore,  in 
view  of  this  rejection  and  the  large  degree  of  skewness,  the  distribution  of 
RMS  values  of  elevation  profiles  filtered  at  500  ft  should  not  be  considered 
as  gaussian. 

3.  1,  2  Power  Spectral  Densities 

The  power  spectral  density  ts  a  functional  characterization 
of  a  time  history  whieh  displays  the  contribution  of  each  wavelength  (or 
frequency)  to  the  total  RMS  of  the  time  history.  Since  the  PSD  not  only 
accounts  for  changes  in  amplitude  of  a  runway  profile  but  also  accounts  for 
the  wavelengths  associated  with  the  amplitude  changes,  such  a  character¬ 
ization  is  important  when  a  runway  is  considered  as  a  dynamic  excitation 
for  an  aircraft.  For  this  reason  the  power  spectral  densities  of  all  of  the 
lines  of  data  were  obtained  and  these  are  presented  in  Appendix  III. 
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Frequency 


Figure  fe  Histogram  of  B-NtS  Value*  With  Profile 
Filtered  at  500  ft 
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Several  comments  are  in  order  regarding  these  figures. 

Power  spectral  densities  had  been  calculated  and  reported  previously  for 
these  data  (c.  f.  Reference  1).  In  these  previous  computations,  the  effect 
of  grades  in  the  profiles  were  removed  by  subtracting  a  linear  mean  as 
defined  by  a  straight  line  between  the  first  and  the  last  points  of  the  profile 
and  the  Fourier  transform  of  the  autocovariance  function  method  was  used 
to  calculate  the  PSD.  Since  such  a  procedure  could  lead  to  a  distortion 
in  the  long  wavelength  region,  all  PSD's  were  rec  'Iculated  for  the  profiles 
digitally  high-pass  filtered  at  X  =  500  and  also  using  the  fast  Fourier 
transform  method  of  computation  as  described  in  Paragraph  2.  3.  In 
general,  the  PSD's  of  this  program  are  in  agreement  with  those  formerly 
calculated  except  for  the  expected  differences  at  the  longest  wavelengths. 

The  differences  for  wavelengths  less  than  200  ft  are  considered  negligible 
in  view  of  the  problem  of  finding  a  general  model  for  all  profiles  at  all 
wavelengths. 

The  PSD's  for  Torrejon  are  distorted  for  the  small  wavelengths. 
This  is  due  to  the  fact  that  the  profiles  for  this  runway  were  taken  by  rod 
and  transit  at  two-foot  intervals  rather  than  the  six-inch  intervals  common 
to  all  of  the  other  profiles.  In  order  to  use  the  same  program  to  process 
all  lines  of  data,  intermediate  points  were  calculated  by  linear  interpolation. 
The  resulting  PSD  is,  of  course,  distorted  for  wavelengths  less  than  four 
feet  with  possible  distortion  due  to  folding  at  higher  wavelengths.  It  should 
also  be  noted  that  the  RMS  values  listed  on  the  figures  of  Appendix  III  are 
those  obtained  as  the  area  under  the  power  spectral  density  as  opposed  to 
those  from  the  time  histories  that  were  discussed  in  the  previous 
paragraph.  For  all  practical  purposes,  the  two  sets  of  RMS  values  are 
in  agreement. 

The  PSD's  of  the  profiles  display  a  reasonable  commonality 
in  shape.  In  the  mid  ranges  of  wavelengths,  they  are  approximately  linear 
(on  a  log-log  scale)  and  they  ail  flatten  at  the  lowest  wavelengths.  The 
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behavior  at  the  longest  wavelengths  is,  however,  erratic.  Many  of  the 
profiles  display  relative  maxima  in  the  mid  ranges  and  most  have  spikes 
at  wavelengths  less  than  five  feet. 

The  flattening  of  the  PSD's  has  generally  been  ignored  in 
modeling  runway  profiles.  It  has  been  conjectured  that  this  flattening  is 
a  result  of  aliasing  from  wavelengths  less  than  one  foot.  To  prove  this 
hypothesis,  a  higher  sampling  rate  than  two  measurements  per  foot  would 
be  required  and  such  data  are  not  available.  Current  specifications  use  a 
linear  extrapolation  from  longer  wavelengths  to  describe  this  region. 

To  investigate  the  gross  commonality  in  shape  of  the  spectral 
densities,  the  PSD  values  at  18  approximately  equally  spaced  reduced 
frequencies  were  normalized  by  the  respective  mean  square  amplitude 
value  for  each  of  the  40  profiles.  This  approach  assumes  that  the  PSD  is 
ol  the  form 


4(0)  =  cr2  $*(0)  (4) 

where  $*(£2)  is  independent  of  a-2  and  is  indicative  of  the  common  shape  of 
all  spectral  densities.  The  resulting  18  normalized  values  for  each  of 
the  40  lines  of  data  were  analyzed  by  determining  the  average,  25th,  ?5th, 
minimum,  and  maximum  values  at  each  of  the  18  reduced  frequencies. 

These  values  are  summarized  in  Figure  7,  For  reduced  wavelengths  less 
than  0.  1  rad/ft,  50  percent  of  the  normalized  spectra  are  within  a  factor 
of  2.  When  engineering  judgment  was  used  to  fit  a  straight  line  (on  log -log 
paper)  to  the  average  normalized  spectra,  the  equation  of  the  line  was 
given  by 

4*(8)  s  0.0125  ST2*00  (5) 

The  scatter  of  the  normalized  spectral  values  is  caused  by 
the  inability  of  Equation  (4)  to  model  the  spectral  densities  and  at  least 
part  of  this  inadequacy  is  the  variation  in  the  slopes  of  the  spectral 
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densities  between  profiles.  To  further  investigate  the  slope  variation,  a 
straight  line  was  fitted  to  the  mid  ranges  of  each  PSD  and  the  slope  was 
determined  for  each  of  the  fitted  lines.  Duo  to  the  engineering  judgment 
involved,  such  slopes  should  bo  considered  accurate  to  within  ±10  percent. 
The  histogram  of  the  slope  values  thus  obtained  are  presented  in  Figure  8. 
The  mean  and  standard  deviation  of  these  values  are  2.  13  and  0,  27, 
respectively.  Thus,  although  the  slopes  do  tend  to  cluster  about  a  value 
of  2,  there  is  significant  variation  in  this  parameter  between  runways. 

It  is  interesting  to  note  the  agreement  between  the  average  slopes  of  2.  13 
and  the  slope  calculated  from  the  mean  of  the  standardized  values  of  2,00 
(c.  f.  Equation  (5)). 

The  erratic  behavior  of  the  spectral  densities  in  the  long 
wavelength  region  can  be  attributed  to  the  inconsistencies  of  the  presence 
of  long  wavelength  undulations  present  in  some  of  the  runway  profiles. 
These  will  be  discussed  further  in  Paragraph  3.  1.  3. 

To  investigate  the  relative  maxima  apparent  in  the  PSD's, 
each  peak  that  was  outside  the  80  percent  confidence  limits  of  adjacent 
points  was  identified  for  all  profiles.  It  must  be  admitted  there  is  a 
degree  of  subjectivity  involved  in  identifying  peaks  in  this  fashion  but 
the  process  tends  to  miss  peaks  rather  than  include  those  which  are 
insignificant.  The  peaks  identified  occurred  at  many  distinct  wavelengths 
up  to  100  ft  but  there  were  clusters  of  peaks  in  certain  wavelength  ranges. 
These  are  summarized  in  Table  V.  According  to  Reference  1,  the  servo 
system  characteristics  of  the  profile  instrumentation  system  coupled 
with  a  constant  surveying  speed  produces  spurious  spectral  peaks  in  the 
two  to  three  foot  wavelength  region.  It  is  speculated  that  the  common 
peaks  apparent  at  multiples  of  2  and  2,5  ft  wavelengths  are  due  to  higher 
harmonies  of  these  spurious  peaks.  The  common  appearance  of  peaks  at 
a  wavelength  of  25.  b  ft  is  probably  due  to  the  tar  strips  or  spalling  between 


TABLE  V 

NUMBER  OF  PROFILES  WITH  RELATIVE  MAXIMA  OF  PSD 
AT  INDICATED  RANGES  OF  WAVELENGTH 


Wavelength 

Number  of  Profiles 

51.  2 

4 

25.6 

16 

15.  1 

10 

12.  5 

9 

10.0 

8 

8.4 

9 

7.  5 

7 

5.0 

6 

4.  0  -  4.  2 

32 

2.  4  -  2.  6 

18 

2.  0  -  2.  1 

36 

1.  3  -  1.4 

42 

the  slabs  of  the  pavement.  If  this  is  true,  a  general  power  spectral  rnodel 
should  contain  such  a  peak  in  this  wavelength  region. 

Two  other  methods  of  considering  the  joint  relationship  of 
amplitude  and  wavelength  were  considered.  Even  though  the  results  from 
these  approaches  were  insignificant,  they  are  described  briefly  in  the 
following  paragraphs. 

The  power  spectral  densities  of  the  raw,  unfiltered  profiles 
were  obtained  and  the  anticipated  results  were  apparent.  The  description 
of  a  ramp  function  in  the  frequency  domain  requires  power  at  ail  wavelengths. 
Therefore,  the  long  term  grades  in  the  original  data  distorted  the  spectral 
densities  by  distributing  the  much  larger  RMS  values  over  all  wavelengths 
to  such  an  extent  that  peaks  in  the  wavelengths  of  interest  were  often  lost. 
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Figures  9  and  10  display  the  PSD's  for  Altus,  Centerline,  and  Bakalar, 
Centerline.  As  shown  in  Reference  1,  Altus,  Centerline,  has  a  linear 
grade  the  entire  length  of  the  runway  with  a  total  change  in  amplitude  of 
approximately  350  inches.  The  PSD  of  this  unfiltered  profile  is  essentially 
a  straight  line.  The  PSD  for  Bakalar,  Centerline,  does  show  vax-iation 
about  a  straight  line  and  the  linear  grade  for  this  profile  extends  only 
22  inches  over  the  length  of  the  runway.  It  may  be  of  interest  to  note  that 
all  of  the  power  spectral  densities  of  the  unfiltered  profiles  could  be 
modeled  quite  well  with  a  straight  line  of  -2  slope  on  the  log-log  plots. 

The  possibility  of  modeling  the  filtered  profiles  by  a  Fourier 
analysis  was  also  briefly  considered.  In  this  approach,  a  profile  was 
assumed  to  represent  one  cycle  of  length  2tt  and  Fourijr  coefficients  for 
wavelengths  between  2  and  512  ft  were  obtained.  The  Fourier  coefficients 
thus  obtained  displayed  such  tremendous  variability  (2  to  3  orders  of 
magnitude)  between  runways  that  this  characterization  was  abandoned. 

3.  1.  3  Stationarity 

An  extremely  important  property  in  the  analysis  of  time 
series  is  that  of  stationarity.  This  property  is  a  necessary  condition  for 
obtaining  statistical  descriptions  from  a  single  time  history.  In  essence  a 
time  series  is  considered  stationary  when  the  statistical  characteristics 
do  not  change  with  time.  Strictly  speaking,  this  property  can  only  be 
established  by  analysis  of  a  collection  of  time  histories  but  in  practice 
this  is  often  impossible.  In  the  analysis  of  the  runway  elevation  profiles, 
it  is  not  sensible  to  attempt  to  establish  stationarity  by  considering  the 
40  lines  of  data  to  be  a  representative  sample  of  the  ensemble  since  known 
differences  exist  in  the  profiles.  For  example,  they  have  different  RMS 
values.  Therefore,  the  investigation  of  stationarity  must  be  performed  on 
the  individual  profiles.  If  the  statistical  properties  do  not  change  with 
distance  down  the  runway,  then  the  profile  c«m  be  said  to  be  self¬ 
stationary  which  is  a  necessary  but  not  sufficient  condition  for  stationarity. 
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REDUCED  FREQUENCY  (R«0/FT) 

Figure  9  Power  Speetral  Deasity  of  UafiUered  Profile  from 
AUus,  Ceaterliue 


RIQUCCO  FRfQUFNCY  (RPQ/FTJ 

Figure  10  Paws?  Spectral  Density  of  Unfittereti  Profile  from 
Bakalar,  Centerline 


Ill  practice,  if  a  time  history  is  shown  to  be  self- stationary,  it  is  treated  as 
though  it  were  stationary  since  quite  often  stationarity  cannot  be  demonstrated 
and  it  is  an  extremely  useful  property. 

The  initial  study  of  stationarity  involved  calculating  the  PSD 
of  each  third  of  each  runway.  If  the  profiles  were  self- stationary,  the 
three  PSD's  for  each  would  be  approximately  the  same.  Part  of  the 
motivation  for  this  approach  was  the  inability  of  the  PSD  method  to  yield 
an  acceptable  prediction  of  a  response  parameter.  One  hypothesis  for  the 
failure  of  the  statistical  approach  to  the  analysis  of  response  was  that  if 
different  segments  of  a  runway  have  different  PSD’s,  then  a  PSD  obtained 
from  the  total  runway  would  not  be  representative  of  the  individual  segments. 

For  several  of  the  runways,  the  three  segments  did  have 
different  PSD's.  As  the  extreme  example,  Figures  11,  12  and  13  present 
the  PSD's  for  the  three  segments  of  the  McGuiro,  Centerline  profile. 

The  decrease  in  power  at  the  long  wavelengths  is  apparent  even  though 
there  is  general  agreement  in  the  shorter  wavelengths.  Note  also  the 
changing  RMS  values  which  decrease  significantlv  between  the  first,  second 
and  third  segments  of  the  runway. 

The  changing  intensity  of  the  remaining  runways  are 
summarized  in  Table  VI  which  displays  the  mean  square  value  for  each 
third  and  the  composite  of  all  three  for  most  of  the  profiles.  Mean  square 
values  are  presented  rather  than  RMS  values  since  they  more  realistically 
reflect  changes  in  the  power  spectral  densities.  It  can  be  seen  that  segment 
mean  square  values  can  disagree  with  the  composite  mean  squares  by  a 
factor  of  8  and  ratios  of  2  to  1  are  not  uncommon.  This  observation  led  to 
a  more  detailed  study  of  the  RMS  values. 


35 


RMS  s  0.  8692 


Figure  12  Power  Spectral  Density  of  Middle  Third  of  MeQuire.  Centerline 
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MEAN  SQUARE  VALUES  FOR  SEGMENTS  OF 


PROFILES  FILTERED  AT  500  FT 


Profile 

Identification 

First 

Third 

Middle 

Third 

Last 

Third 

Composite 

Altus 

75  L 

mm 

0.  109 

0.  168 

0.  131 

CL 

0.014 

0.  0S3 

0.033 

Bakalar 

50L 

0.  150 

0.  196 

0.  209 

0.  195 

CL 

0.  060 

0.  066 

0.079 

0.070 

Carswell 

37L 

0.  114 

0.033 

0.087 

13L 

0.  130 

0.081. 

0.028 

0.078 

13R 

0.  153 

0.  096 

0.039 

0.095 

Edwards 

75L 

0.  087 

0.047 

0.044 

0.060 

CL 

0.041 

0.045 

0.045 

0.044 

75R 

0.  054 

0.  038 

0.036 

0.043 

Ellsworth 

CL 

0.039 

0.073 

0.  122 

0.076 

Glasgow 

75L 

0.077 

0.  142 

0.  107 

0.  109 

CL 

0.045 

0.  106 

0.065 

0.072 

Largley 

SOL 

0.077 

0.  165 

0.  172 

0.  136 

CL 

0,070 

0.  123 

0.084 

0.089 

McGuire 

75L 

1.  168 

0.  543 

0.  140 

0.620 

CL 

2,648 

0.756 

0.  149 

1.  201 

75R 

0.  574 

0.610 

0.  202 

0.454 

Shaw 

SOL 

0.031 

0.060 

0.  105 

0.064 

Tr-avis 

75L 

0.415 

0.  243 

0.470 

0.  361 

CL 

0.  183 

0.  134 

0.093 

0.  135 

75R 

0.599 

0.  364 

0.  258 

0.430 

Torrejon 

12L 

0.  201 

0.084 

0.  032 

0.  106 

CL 

0.  236 

0.093 

0.069 

0.  132. 

12R 

0.  175 

0.082 

0.047 

0.  101 

Castle 

75L 

0.  345 

0.  196 

0.  134 

0.  220 

CL 

0.048 

0.087 

0.077 

0.069. 

75R 

0.  138 

0.403 

0.  110 

0.  220 

George 

50L 

0.  121 

0.072 

0.  105 

0.097 

CL 

0.  115 

0.  101 

0.  133 

0.  115 

Palmdale 

50L 

0.  102 

0.084 

0.072 

0.087 

CL 

0.068 

0.  108 

0.070 

0.082 

50R 

0,064 

0.094 

0.078 

0.077 

The  RMS  values  for  each  512  ft  section  of  each  profile  were 
obtained  and  plotted  as  a  function  of  distance  down  the  runways.  These  RMS 
histories  for  all  of  the  profiles  are  presented  in  Appendix  IV  and,  as  an 
example,  the  RMS  histories  for  Altus  and  Bakalar  are  repeated  as  Figure  14. 
It  is  apparent  that  although  some  of  the  profiles  display  a  consistency  of 
RMS  values  over  the  entire  length  of  the  runway,  many  also  display 
significant  changes  between  segments.  Therefore,  it  is  concluded  that  the 
assumption  of  stationarity  is  not  true  in  general  even  though  some  of  the 
profiles  are  self-stationary. 

To  further  investigate  the  high  RMS  regions  of  the  profiles, 
the  corresponding  regions  were  identified  on  plots  of  the  runway  profiles. 

In  every  case,  the  high  RMS  values  could  be  attributed  to  high  amplitude, 
long  wavelengths  undulations  as  opposed  to  intense  short  wavelength  regions 
of  roughness.  This  observation  leads  naturally  to  the  hypothesis  that  the 
long  wavelength  undulations  are  the  cause  of  the  nonstationarity. 

3.  2  MAXIMUM  WAVELENGTH  OF  100  FEET 

Since  the  analysis  of  the  profiles  filtered  at  a  maximum  wavelength  of 
500  ft  indicated  that  the  long  wavelength  undulations  are  the  possible  cause 
of  both  nongaussianness  and  nonstationarity,  the  decision  was  reached  to 
high-pass  filter  all  profiles  at  a  lower  maximum  wavelength  and  to  analyze 
these  different  elevation  profiles.  Of  course,  the  longer  wavelengths  are 
important  ar  dynamic  excitation  to  aircraft  and  they  must  be  accounted  for 
in  the  characterization  of  runway  profiles.  However,  in  the  synthesis  of  a 
simulated  runway,  it  would  be  most  convenient  if  the  shorter  wavelengths 
could  be  modeled  in  the  frequency  domain.  Since  the  previous  analysis  has 
shown  that  all  wavelengths  up  to  500  ft  cannot  be  modeled  in  the  frequency 
domain,  a  compound  characterization  of  runway  roughness  would  be 
required.  The  long  wavelengths  would  be  modeled  by  a  discrete  description 
of  amplitudes  and  wavelengths  whereas  the  shorter  wavelengths  would  be 


modeled  by  a  gaussian,  stationary  time  history  with  a  PSD  as  determined 
by  actual  runway  profiles. 

Two  considerations  entered  into  the  choice  of  the  reduced  maximum 
cutoff  wavelength.  For  ease  in  modeling,  it  is  desirable  to  have  as  large 
a  maximum  value  as  possible  for  which  stationarity  and  gaussianness  can 
be  reasonably  assumed.  On  the  other  hand,  since  the  dominant  excitation 
occurs  at  twice  the  gear  spacing  of  the  aircraft,  the  choice  of  the  maximum 
wavelength  for  the  short  wavelength  characterization  should  be  such  that 
it  would  contain  wavelengths  of  twice  the  gear  spacings  of  large  aircraft. 

For  these  reasons,  the  initial  choice  of  a  reduced  maximum  wavelength  was 
chosen  at  100  ft  with  the  intention  of  increasing  this  maximum  if  possible. 

As  will  be  shown,  however,  stationarity  and  gaussianness  are  only 
marginally  acceptable  for  all  wavelengths  less  than  100  ft  so  that  a  larger 
maximum  cutoff  would  lead  to  a  higher  degree  of  nonstationarity  and  non- 
gaussianness.  Therefore,  no  other  maximum  wavelength  was  considered. 

All  profiles  were  digitally  high- pass  filtored  by  the  methods  described 
in  Paragraph  2.  2.  All  wavelengths  less  than  100  ft  were  passed  unattenuated, 
all  wavelengths  greater  than  140  ft  were  essentially  eliminated,  and  between 
100  and  140  ft  the  amplitudes  were  attenuated  in  accordance  with  the  filter 
response  of  Figure  2.  The  analyses  of  Pai’agraph  3.  1  were  then  repeated 
for  these  new  profiles  and  the  results  of  these  analyses  are  presented  in 
the  following  paragraphs. 

3.  2.  1  Amplitude  Probability  Densities 

The  amplitude  probability  densities  were  obtained  for  each  line 
of  data  filtered  at  a  100  ft  maximum  wavelength  and,  again,  the  cumulative 
distribution  functions  were  plotted  on  gaussian  probability  paper.  These 
comparisons  for  ail  lines  of  data  are  presented  in  Appendix  V.  Again,  the 
solid  line  represents  the  gaussian  fit  as  obtained  from  the  mean  and  standard 
deviation  of  the  amplitudes  and  not  a  best  engineering  judgment  line  through 
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the  points.  For  all  lines,  the  amplitude  probability  densities  were  more 
gaussian  than  those  obtained  from  filtering  at  a  max' mum  cutoff  vavelength 
of  500  ft.  However,  discrepancies  in  the  tails  of  the  distributions  are  still 
present  though  less  pronounced  than  before.  Some  of  the  extreme 
discrepancies  may  be  explainable  by  noise  spikes  which  passed  the  editing 
criteria  of  0.4  inches  described  in  Paragraph  2.  1.  The  RMS  values  of  the 
lines  filtered  at  100  ft  are  considerably  smaller  (average  value  of  0.  16  in) 
than  those  filtered  at  500  ft  (avex'age  value  of  0.  36  in).  Therefore,  the 
effect  of  these  spikes  may  have  a  greater  effect  in  the  distributions  of 
amplitudes  filtered  at  100  ft.  Again,  the  gaussian  distribution  would  provide 
the  best  single  family  of  distributions  in  modeling  the  amplitude  probability 
density. 

The  RMS  values  of  the  lines  filtered  at  100  ft  are  presented 
in  Table  VII  for  each  profile.  These  values  are,  of  course,  lower  than  the 
corresponding  RMS  value  from  the  lines  filtered  at  500  ft.  A  histogram  of 
the  reduced  RMS  values  is  presented  in  Figure  15  from  which  it  is  easily 
seen  that  the  scatter  in  the  RMS  values  is  also  considerably  reduced  as 
compared  to  tho>s  from  the  profiles  filtered  at  500  ft.  For  these  RMS 
values  the  mean,  standard  deviation,  skewness  and  kurtosis  are  given  by 
0.156,  0.042,  1.  992  and  6.  916,  respectively,  A  Kolmogoroff-Smirnov 
goodness  of  fit  test  rejects  the  hypothesis  that  these  RMS  values  are  from 
a  normal  distribution  at  the  99  percent  level  of  confidence.  It  is  not 
necessary  for  these  RMS  values  to  be  normally  distributed  but  if  design 
criteria  for  intensity  are  set  by  the  mean  plus  three  standard  deviations 
of  the  RMS  distribution,  there  is  no  assurance  that  this  value  will  be  greater 
than  99.  9  percent  of  the  intensities  to  be  ineouotered.  Actually,  the  cause 
of  the  rejection  of  normality  is  the  three  large  RMS  values  between  0.  <i75 
and  0.  300  which  came  front  the  lines  of  data  for  McGuire, 
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3.  2.  2  Power  Spectral  Densities 


Theoretically,  there  should  be  no  difference  between  the  power 
spectral  densities  of  profiles  filtered  at  100  ft  and  those  filtered  at  500  ft  in 
the  range  of  wavelengths  less  th.v.n  100  ft.  To  check  this  theoretical  result, 
the  PSD's  were  recomputed  for  -ill  of  the  lines  with  the  reduced  cutoff 
wavelength.  Since  these  PSD1?  agreed  very  closely  for  wavelengths  less 
than  100  ft  with  those  already  presented  in  Appendix  III,  all  of  the  comments 
of  Paragraph  3.  1.  2  regarding  analysis  of  the  PSD's  in  the  mid  ranges  of 
wavelength  are  directly  applicable.  Therefore,  no  further  analysis  is 
required.  Perhaps  it  sho  tld  be  noted  that  the  cutoff  wavelength  of  100  ft 
eliminates  only  the  first  four  points  of  the  spectra  of  Appendix  III  with  a 
slight  attenuation  of  the  fifth  point  which  is  at  a  wavelength  of  102  ft. 

3.  2.  3  Station?  rity 

To  in’  estigate  stationarity  of  the  profiles  filtered  at  100  ft, 

RMS  values  were  again  calculated  for  each  512  ft  section  of  each  profile 
and  plotted  as  a  function  of  distance  down  the  runway.  These  data  are 
presented  in  Appendix  VI  and,  for  illustrative  purposes,  the  RMS  histories 
for  Altus  an  i  Bakalar  are  repeated  in  Figure  16.  It  is  apparent  from  these 
figures  th'.t  the  degree  of  self-stationarity,  as  indice.ted  by  RMS  values,  is 
much  be. ter  for  the  profiles  filtered  at  100  ft  than  at  500  ft.  Approximately 
half  of  the  lines  display  erratic  behavior  in  RMS  and  these  generally  have 
isolated  sections  of  one  RMS  value  considerably  higher  than  the  others. 

In  most  of  the  lines,  however,  the  erratic  behavior  has  been  significantly 
.•educed  by  the  filtering  at  a  lower  cutoff  wavelength.  The  three  lines  from 
McGuire  again  standout  as  being  different  from  the  remaining  lines. 

In  the  compound  characterization  of  runway  roughness,  the 
degree  of  non-self-stationarity  would  have  less  effect  in  the  ayntheais  of  a 
simulated  runway.  This  lessening  of  effect  would  be  due  to  the  reduced 
RMS  values  for  the  short  wavelengths  coupled  with  the  independent 
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introduction  of  the  long  wavelength,  high  amplitude  cycles.  Further,  it  was 
not  considered  practical  to  go  to  a  maximum  wavelength  shorter  than  100  ft. 
Therefore,  even  though  there  is  evidence  of  nonstationarity,  it  was  assumed 
in  modeling  the  wavelengths  less  than  100  ft,  that  the  short  wavelength 
contribution  to  the  simulated  runways  could  be  synthesized  by  using  a 
theory  dependent  on  the  assumption  of  stationarity. 

3.  3  DISCRETE  LONG  WAVELENGTH  DESCRIPTION 

If  a  compound  characterization  of  runway  profiles  is  adopted  wherein 
the  wavelengths  less  than  100  ft  are  modeled  in  the  frequency  domain  using 
available  thoory  for  stationary,  gaussian  time  histories,  the  longer  wave¬ 
length  contribution  to  runway  profiles  must  be  independently  characterized. 
Two  methods  of  characterization  of  the  long  wavelengths  components  were 
giv  rn  serious  consideration.  Initially,  relative  maxima  and  minima  were 
paired  by  an  algorithm  with  the  distance  between  these  peaks  defined  as  the 
half  wavelength  and  the  difference  in  magnitudes  defined  as  twice  the 
amplitude.  When  the  algorithm  was  applied  to  the  profiles  filtered  at 
500  ft,  however,  many  of  the  peaks  were  described  as  having  wavelengths 
greater  than  500  it.  Since  such  a  result  is  contradictory,  this  method  of 
characterization  was  abandoned. 

Study  of  the  individual  histories  of  the  profiles,  indicated  that  the 
long  wavelength  components  are  not  generally  superimposed.  That  is, 
when  secondary  peaks  are  present  on  a  long  wavelength  undulation,  these 
are  due  to  cycles  of  less  than  100  ft  wavelength.  Therefore,  the  method 
adopted  far  characterizing  the  long  wavelength  undulations  was  that  of 
primary  peaks  of  the  profiles.  Each  maximum  and  each  minimum  between 
zero  crossings  was  identified  and  associated  with  a  half  wavelength  defined 
as  the  distance  between  aero  crossings.  All  peaks  with  a  half  wavelength 
less  than  SO  ft  were  then  dropped  as  these  will  be  modeled  by  the  short 
wavelengths  component  of  runway  profiles. 
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The  characterization  of  the  long  wavelength  undulations  is  then 
expressed  in  terms  of  the  distribution  of  peaks  per  1000  ft,  the  distribution 
of  half  wavelengths,  and  the  joint  distribution  of  half  wavelengths  and 
amplitudes.  The  method  for  using  these  distributions  in  deriving  a 
simulated  runway  is  presented  in  the  next  section. 

Since  the  available  runway  profiles  were  of  different  lengths,  and 
since  it  may  be  desirable  to  simulate  runways  of  different  lengths,  the 
number  of  primary  peaks  in  each  profile  was  normalized  by  the  length  of 
the  runway.  The  cumulative  distribution  of  peaks  per  1000  ft  was  then 
obtained  and  this  distribution  is  presented  in  Figure  17  on  normal 
probability  paper.  Although  the  hypothesis  of  normality  could  not  be 
rejected  by  the  Kolmogoroff-Smirnov  test,  the  actual  observed  distribution 
was  used  in  the  runway  simulation  process  rather  than  a  normal  distribu¬ 
tion  fit  of  data  points.  For  the  distribution  of  peaks  per  1000  ft  the  mean 
and  standard  deviation  are  given  by  5,  26  and  0.  69,  respectively. 

In  this  simulation  process  knowing  the  number  of  peaks  per  1000  ft 
is  not  sufficient  without  knowledge  of  the  manner  in  which  these  are 
distributed  along  the  runway.  As  each  peak  was  identified,  its  location  on 
the  runway  was  also  noted.  For  each  individual  line  of  data,  the  number 
of  peaks  in  1000  ft  increments  down  the  runway  was  obtained.  Each  of  these 
40  distributions  were  tested  against  a  uniform  distribution  for  location  of 
the  peaks  by  a  Kolmogoroff-Smirnov  goodness  of  fit  test  at  a  level  of 
confidence  of  80  percent.  The  equality  of  the  observed  distribution  of 
location  of  peaks  to  the  uniform  distribution  was  rejected  in  only  tnree  of 
the  40  tests.  Since  testing  was  performed  at  a  level  ot  confidence  of 
80  percent,  the  expected  number  of  rejections  in  40  tests  is  eight  if  the 
distribution  of  peak  location  is  truly  uniform.  Therefore,  these  tests 
indicate  that  the  peaks  are  uniformly  distributed  over  the  lengths  of  the 
runways. 


The  cumulative  distribution  of  half  wavelengths  is  presented  in 
Figure  18  and  the  cumulative  distributions  of  amplitudes  for  the  half 
wavelength  intervals  are  presented  in  Table  VIII.  In  the  generation  of  the 
simulated  runway  these  observed  distributions  wore  directly  used  rather 
than  fitting  them  with  some  family.  The  dependence  of  the  amplitude 
distributions  on  half  wavelength  is  apparent  from  Table  V5II.  Note,  for 
example,  the  increasing  median  value  with  wavelength.  Nofe  also  the 
increased  proportion  of  high  amplitude  values  at  the  long  wavelength 
intervals.  This  dependence  is  accounted  for  in  the  simulation  process 
by  first  drawing  a  wavelength  at  random  from  the  cumulative  distribution 
of  Figure  18  and  then  drawing  an  amplitude  at  random  from  the  appropriate 
cumulative  distribution  of  Table  VIII.  This  procedure  is  an  application  of 
the  well* known  law  of  probability  that 

P  {A  and  B  }  P  { B  }  Pf  A  given  B } 


where 

A  denotes  a  primary  peak  of  given  amplitude 
B  denotes  a  primary  peak  of  given  wavelength 

This  procedure  insures  a  random  selection  of  peaks  with  representative 
amplitudes  and  wavelengths. 
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SECTION  4 

SIMULATED  RUNWAYS 


A  principal  result  of  the  analysis  of  the  runway  roughness  data  was 
to  produce  a  simulated  runway  for  use  as  design  criteria.  This  simulated 
runway  must  have  statistical  properties  similar  to  the  sample  of  discrete 
runway  data.  A  decision  was  made,  during  the  analysis  of  the  data,  to 
represent  the  relatively  high  frequency  (less  than  100  ft  wavelength) 
component  in  terms  of  a  PSD  and  the  low  frequency  data  in  terms  of  discrete 
peaks  and  valleys.  Using  the  following  equation  {Reference  2)  a  time 
history  can  be  computed  from  any  PSD. 


where 


N 

f(x)  =  Yj  c os(x  ujj  -  0j) 

i=l 

f(x)  is  the  amplitude  at  point  x  on  the  runway 
At  =  (c^OD)  AU))1/2 

<j>^((D)  is  the  PSD  as  a  function  of  frequency 
is  the  frequency 

8i  is  a  random  angle  between  0  and  2n 
N  is  the  number  of  points  used  to  describe  <jq(U)) 


(6) 


Using  Equation  (6)  a  runway  of  any  length  can  be  created  and  it  will  have 
the  desired  PSD,  hence  the  same  RMS  value.  Large  wavelength  peaks  and 
valleys  can  then  be  superimposed  on  this  profile  as  desired. 
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4.  1  SAMPLE  RESULT 


A  computer  program  has  been  written  to  generate,  at  random, 
typical  runways  with  the  high  frequency  content  described  as: 


§(ft)  = 


.0023 


a 


2.007 


(7) 


where 


$(G)  is  the  PSD  (in2/(rad/ft)) 

£2  is  the  reduced  frequency  (rad/ft) 

This  equation  was  derived  by  fitting  a  straight  line  (on  log-log  paper)  to 
the  mean  PSD  of  all  lines  of  survey  to  obtain  the  slope  and  using  the  RMS 
value  for  which  90  percent  of  the  lines  of  survey  do  not  exceed. 

The  low  frequency  content  was  added  to  the  runway  by  first  randomly 
selecting  the  number  of  peaks  to  be  added  from  Figure  17,  which  is  the 
distribution  of  long  wavelength  peaks  obtained  from  the  runway  profiles. 
This  distribution  contains  the  sum  of  the  positive  and  negative  peaks  since 
it  was  found  that  there  were  an  equal  number.  The  rate  drawn  from  the 
distribution  is  then  multiplied  by  10  for  the  10,  000  ft  runway.  Next,  for 
each  peak,  which  are  alternately  assigned  positive  and  negative  values,  a 
half  wavelength  was  randomly  drawn  from  the  cumulative  distribution  of 
half  wavelengths  (Figure  18)  and  an  amplitude  was  randomly  drawn  from 
the  cumulative  distribution  of  amplitudes  corresponding  to  the  wavelength 
(Table  VIII).  The  long  wavelength  peaks  were  determined  to  be  uniformly 
distributed  over  the  length  of  the  runway,  however,  selecting  a  random 
place  for  the  peaks  becomes  a  matter  for  discussion.  After  some 
deliberation,  it  was  decided  that  the  peaks  would  be  superimposed,  as 
1/2  sine  waves,  at  the  center  of  a  randomly  selected  available  segment  of 
the  runway.  Therefore,  the  first  peak  selected  is  put  in  the  center  of  the 
runway,  leaving  two,  equal  length  available  segments.  The  second  peak 
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is  then  put  in  the  center  of  one  of  the  two  available  segments,  which  now 
leaves  three  available  segments.  This  procedure  is  continued  until  all 
peaks  have  been  superimposed  on  the  runway.  If  at  any  time  the  randomly 
selected  segment  is  too  short  to  accept  the  peak,  successive  segments  are 
examined  until  a  long  enough  one  is  found,  or  if  none  are  long  enough  a 
message  is  written  by  the  program  and  the  peak  is  not  included. 

Interpolation  was  performed  in  determination  of  the  wavelength  and 
amplitude  of  each  peak,  that  is,  the  distributions  were  considered  to  be 
continuous  from  50  ft  to  300  ft  in  wavelength  and  from  0.  25  inch  to  5.  0 
inches  in  amplitude.  The  number  of  peaks  per  thousand  feet  of  runway, 
however,  was  not  interpolated  since  there  are  36  entries  in  the  table  and 
the  span  when  multiplied  by  10  is  only  35. 

In  order  to  prove  that  a  runway  generated  in  such  a  manner  has  the 
desired  statistical  properties,  the  simulated  runway  was  then  input  to  the 
program  which  extracts  the  primary  peaks  from  the  data.  TabLes  IX  and 
X  show  the  peaks  that  were  superimposed  on  the  runway  and  the  peak 
distributions  that  were  extracted  from  the  runway,  respectively.  The 
profile  was  then  filtered  using  a  100  ft  cutoff  and  the  PSD  was  computed, 

A  plot  of  this  PSD  is  presented  in  Figure  19  and  the  corresponding 
amplitude  probability  density  in  Figure  20.  The  entire  profile  of  the 
simulated  runway  is  presented  in  graphical  form  in  Figure  21  and  in 
tabular  form  in  Table  XI. 
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TABLE  IX 

DISTRIBUTION  OF  LONG  WAVELENGTH  PEAKS 
INPUT  TO  SIMULATED  RUNWAY 


Amplitude 


0.  375 
0.  500 
0.625 
0.  750 
0.875 
1.000 
1.  125 
1.  250 

1.  375 
1.500 
1.625 
1.750 
1.875 

2.  000 

2.  250 


j  Wavelength 

70 

90 

110 

130 

150 

200 

250 

6 

2 

2 

6 

4 
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6 

1 

2 

1 

5 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

I 

1 

1 

1 

1 

1 

TABLE  X 

DISTRIBUTION  OF  LONG  WAVELENGTH  PEAKS 
OUTPUT  FROM  SIMULATED  RUNWAY 


Amplitude 


.  375 
.  500 
.625 
.750 
0.  875 
1.000 
1.  125 
l.  250 

1.  375 
1.500 
1.625 
1.750 
1.875 

2.  000 

2.  250 


Wavelength 

70 

90  110  130  150 

200 

250 

1  1 

2  1 
1  1 
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Figure  19  PSD  from  Simulated  Runway 
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4.  2  OTHER  USES  OF  SIMULATION  PROGRAM 


The  program  written  to  generate  simulated  runway  profiles  can  be 
altered  to  perform  simulations  of  many  different  types  of  runways.  Any 
PSD  may  be  input  to  that  portion  of  the  program  which  generates  data  from 
the  PSD.  T.'ae  sample  result  uses  a  linear  PSD  to  generate  data  with 
wavelengths  of  100  ft  down  to  one  foot,  however,  this  could  easily  be 
changed  to:  1)  cover  a  different  band  of  wavelengths,  2)  represent  a  PSD 
with  a  different  slope  or  RMS,  or  3)  input  any  given  PSD  from  cards.  The 
method  of  generating  discrete  long  wavelength  peaks  could  also  be  alterod 
to:  1)  generate  a  givon  number  of  peaks,  2)  generate  a  profile  containing 
specific  peaks  at  specific  places  on  the  profile,  or  3)  eliminate  the 
generation  of  long  wavelength  peaks  completely. 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


Runway  elevation  profiles  of  40  lines  of  data  from  >4  runways  were 
analyzed  to  achieve  two  objectives:  1)  to  determine  and  evaluate  significant 
characteristics  of  runway  profiles  with  respect  to  their  influence  as  dynamic 
excitation  of  aircraft  and  2)  to  develop  a  method  for  generating  simulated 
runways  which  have  these  significant  characteristics.  To  achieve  these 
joint  objectives,  runway  elevation  profiles  of  40  lines  of  data  from  14 
runways  were  assumed  to  be  a  representative  sample  of  the  population  of 
all  runways  of  interest  and  analyzed  with  both  objectives  in  view. 

The  grades  present  in  the  raw  data  were  first  removed  by  a  digital 
high- pass  filter  with  a  maximum  cutoff  wavelength  of  500  ft.  The 
resulting  profiles  were  shown  to  be  both  nunstationary  and  nongaussian. 

The  failure  of  the  profiles  to  meet  these  basic  assumptions  provides  one 
explanation  for  the  failure  of  the  statistical  approach  to  the  prediction  of 
the  response  and  also  precluded  the  possibility  of  generating  a  simulated 
runway  by  modeling  in  the  frequency  domain  over  all  wavelengths  less 
than  500  ft.  Further  analysis  of  the  causes  of  nonstationarity  and  nen- 
gaussianneas  indicated  that  to  a  targe  extent  they  are  the  result  of  the 
larger  wavelength  undulations  of  the  profiles. 

The  data  were  again  high- pass  filtered  with  a  reduced  maximum 
wavelength  of  100  ft  and  for  these  new  profiles  the  degree  of  stationarity 
and  gaussiaaness  was  significantly  better  than  for  the  profiles  filtered  at 
500  ft.  Since  the  shorter  wavelength  profiles  have  a  smaller  intensity 
(RMS)  than  the  longer  wavelength  profiles,  the  effects  of  neostatiaaarity 
and  nongaussianness  in  a  characterisation  of  a  runway  with  all  wavelengths 
up  to  500  ft  would  be  minor.  Therefore,  a  compound  characterisation  was 
derived  wherein  the  wavelengths  less  than  100  ft  are  modeled  in  frequency 
domain  and  the  longer  wavelengths  are  modeled  in  a  discrete  fashion  by 


?8 


means  of  the  distributions  of  primary  peaks  per  1000  ft  and  the  joint 
distribution  of  wavelengths  and  amplitudes  of  the  primary  peaks.  Finally, 
a  method  is  presented  for  generating  a  simulated  runway  which  has  the 
properties  of  the  compound  characterization  of  the  actual  profiles. 

it  should  be  stressed  that  the  PSU  plot  shown  in  Figuro  19  is  for 
the  filtered,  simulated  runway.  As  can  bo  seon  in  this  plot,  the  region 
groator  than  one  hundred  foot  is  greatly  attenuated.  The  filter  attenuates 
all  amplitudes  with  a  wavelength  greater  than  one  hundred  feet  and  this 
PSD  plot  is  only  valid  for  wavelengths  less  than  one  hundred  feet. 

The  runway  profile  elevations  presented  in  Table  XI  merely 
demonstrates  the  use  of  ih.  technique  developed  in  this  study  and  they 
are  not  intended  for  aircraft  design  purposes. 


T9 


APPENDIX  I 

DIGITAL  FILTERING  CONCEPTS 


The  digital  filtering  process  consisted  of  subtracting  the  low  frequency 
* 

mean,  obtained  by  applying  a  low-pass  filter  to  the  original  data,  from  the 
original  data.  This  yields  data  which  lias  been  high-pass  filtered  and  has  a 
zero  mean.  The  equation  illustrating  the  high-pass  filter  operation  is 

*d  =  xd  "  x  d  (3) 

where 

x^  is  the  sampled  input 

x'^  is  the  output  of  the  low-pass  filter 

x^j  is  the  final  high-pass  filtered  data 

The  low-pass  filtering  operation  for  equi-spaced  intervals,  Ad,  is  given  as 


N 

xda  l'  hn*(d+nAd)  <8) 

n=-N 


where 

hM  are  the  2N  +  1  filter  weights 
:<(d+nAd)  are  the  input  data  samples 
X  jj  are  the  output  data  samples 
Ad  is  the  distance  between  samples 


The  filter  shape  chosen  is  attributed  to  Martin  and  Graham  who  developed 
it  independently  (Reference*  3  and  4).  The  weights  for  this  filter  are 
computed  from 


dwnAd 


^$i«2ir(ffc  +  Af)nAd  +  siudiff^nAd 

I  -  4(AfnAd)“ 


(9) 


£0 


where 


n  *  1,  ±2,  ...  ±  N 
fc  is  the  cutoff  frequency 
Af  is  the  roil-off  frequency 
and  ^n=0  =  +  Af 


Certain  trends  given  as  general  polynomial  forms  and  considered  as 
desirable  data  components  require  constraints  on  the  basic  weights  in 
order  to  pass  without  error.  If  the  function  does  contain  a  polynomial, 
content,  it  is  not  band  limited  and  should  be  constrained  so  that  the  sample 
values  of  the  polynomial  content  are  passed  without  error.  Therefore, 
to  prevent  trend  type  errors,  constraints  must  be  applied  to  the  weighting 
function  to  satisfy 

iSaifl=oi 

4t’<s  j  f=Q  for  k  a  1 

(10) 


f=0  for  k  -  0 


For  k  ?  0  as  in  the  eu«e  with  the  Martiu-Uraham  filter,  this  is  accomplished 
by  normalising  the  weights  so  that 

N 

Z  hn«l  01) 

na-N 

Therefore,  the  normalised  weights  are  computed  from 


is-N 


(U) 


SI 


This  forces  the  transfer  function  to  have  unity  gain  at  aero  frequency.  It 
should  be  noted  that  h_n  =  hn,  so  that  only  weights  for  positive  n  need  be 
computed. 

Application  of  the  digital  filter  involves  the  selection  of  fc,  Af  and  N. 
Too  small  a  Af  or  too  small  a  value  of  N  will  cause  excessive  overshoot  at 
the  filter  corner  frequencies.  In  extreme  cases,  this  overshoot  can  carry 
through  a  wide  range  of  the  spectrum.  The  sampling  frequency,  f„  =  1/Ad, 
must  be  greater  than  twice  the  highest  significant  frequency  present  in  the 
data  to  prevent  aliasing.  However,  too  large  a  value  of  fs  will  reduce  the 
distance  width  of  tho  filter  weighting  function.  The  only  remedy  is  to 
increase  N,  which  is  usually  impractical,  or  to  decrease  fs.  This  is 
accomplished  by  pro-filtering  the  data  through  a  Low-pass  filter  that  has 


where 


f,2/2  *<fc,  +*£])/<! +  (fc^  +  Af2) 


(13) 


fCj  =  cutoff  frequency  of  the  preliminary  low-pass  filter 
HCy  =  cutoff  frequency  of  the  final  low-pass  filter 

h~t 

Afj  5  the  preUninary  filter  roll-off 

a  the  sampling  frequency  of  the  data  to  the  final 
tow-pass  filter 

Af^  s  the  final  filter  roll-off 


The  requirements  of  the  initial  filter  were  that  the  frequencies  to  be 
passed  by  the  final  low-pass  filter  remain  unchanged  and  that  the  gain  of 
the  filter  transfer  function  be  aero  at  f^  ±-  (fe?  *  Af,?)*  The  sero  gain 
requirement  is  necessary  so  that  there  is  no  significant  frequencies 
present  above  the  Nyquist  frequency  which  would  introduce  aliasing. 


(14) 


=  2Ad2  =  2 

where  is  the  Nyquist  frequency. 

Aliases  that  occur  at  frequencies  in  the  output  of  the  preliminary  filter 
greater  than  fc  +  Af2  are  of  no  concern  since  they  will  be  filtered  out  by 
the  final  low-pass  filter.  The  expression  for  the  double  filter  is 

0  N?  N> 

xd  =  L  ^n2  Zj  ,,c(d+niAdi+ ti^Ad?)  (*) 

n2=-N2  nj=-Nj 

Filtered  high-pass  data  is  obtained  by  linear  interpolation  between  samples 
of  the  low-pass  data  v/ith  sampling  frequency  fs  and  subtracting  from  the 
original  data  with  sampling  frequency  fSj.  Consequently,  to  keep  interpola¬ 
tion  errors  to  a  minimum,  fs^  must  be  chosen  large  enough  to  reconstruct 
the  highest  frequency  passed  by  the  low  .-pass  filter.  The  following  guidelines 
are  used  to  determine  fs,  Li  and  N,  given  fg,  for  a  resultant  overshoot  S  1 
percent  and  an  interpolation  error  5  1  percent. 


fg  £  20(fe  +  Af) 


5  0.015  «f 


N  & 


1.0Sfg 

Li 


OS) 


One  additional  guideline  must  also  he  fulfilled  in  order  to  avoid  a  limit 
condition  in  the  computation  of  the  filter  weights.  This  is  0. SfLiLd  may 
not  be  an  integer  value  «  N.  The  empirical  error  b  uad  for  the  Martin- 
Graham  filter  in  terms  of  N,  Li  and  &d  is  given  as 


In 


A  “  i  iV  4 


4Na(£f£d)2  -  1 


|*U| 


a* 


An  indication  of  just  how  well  the  two  stage  filters  designed  for  filtering 
the  runway  profile  data  fulfilled  the  suggested  guidelines  and  the  estimates 
of  the  error  associated  with  the  filters  can  be  found  in  Table  XU. 

TABLE  XII 
FILTER  GUIDELINES 


Guideline 

500  ft  Cutoff 

100  ft  Cutoff 

Preliminary 

Filter 

Final 

Filter 

Preliminary 

Filter 

Egg 

fs  >  20  (fc  +  Af) 

0.63 

0.04 

0.8 

0.2 

(2) 

(0.0667) 

(2) 

(0.  2) 

Af  0.0l5fs 

O 

• 

o 

0.001 

0.03 

1 

(0.03) 

(0.0005) 

(0.033) 

■ 

N  >  1. 05fg/Af 

70 

140 

64 

,1 

(75) 

(HO) 

Bl  I 

n 

1.  4% 

1.6% 

1.0% 

0.  ?•% 

NOTE:  Values  in  parentheses  are  the  actual  values  of  f$,  Af  and  N  used. 

In  the  actual  application  of  the  low-pass  filters,  two  different 
computer  programs  were  written.  The  first,  for  the  500  ft  cutoff,  applied 
the  preliminary  filter  in  the  distance  domain.  This  makes  a  quite  complex 
program  in  which  the  filter  constants  cannot  be  changed  by  much  without 
restructuring  the  program.  The  final  stage  was  then  applied  in  the 
frequency  domain  by  making  use  of  the  Fast  Fourier  Transform  and  the 
convolution  theorem  which  states  that 

Y(s)  «  Wftt)X(a)  (17) 

The  interpretation  of  this  expression  is  that  the  Fourier  transform  of  the 
faltered  output  is  equal  to  the  prvduct  of  the  Fsurter  trass  farms  at  the 
filter  transfer  function  and  the  input  function.  Therefore,  the  second  stage 


weights  wore  computed  using  Equations  (9)  and  (12)  and  then  transformed. 

The  output  of  the  first  stage  filter  was  then  transformed,  the  complex  product 
of  the  two  transforms  was  formed  and  the  inverse  transform  was  made 
resulting  in  the  final  low-pass  filtered  data. 

The  advantage  of  applying  the  first  stage  filter  in  the  distance  domain 
is  that  it  takes  less  computer  time  since  it  is  only  necessary  to  compute 
the  output  lr,c  every  30th  point. 

When  it  was  decided  that  the  data  should  be  filtered  with  a  shorter 
wavelength  cutoff,  a  new  filter  program  was  written  to  make  use  of  the  FFT 
in  both  stages  of  the  filter,  I3y  doing  this,  a  more  general  filter  program, 
in  which  the  filter  constants  could  easily  b  *  varied,  was  obtained.  The 
time  increase  of  approximately  15  percent  for  this  program  was  of  little 
significance. 

The  only  additional  aspect  of  the  filtering  procedure  not  yet  discussed 
is  that  of  folding  the  data  at  the  beginning  and  end  of  the  runway  profile. 

When  data  are  filtered,  N  values  are  lost  from  the  beginning  and  end.  of  the 
series.  Several  schemes  are  available  to  extend  the  data  so  that  the  output 
is  the  same  length  as  the  input.  The  choice  of  scheme  to  use  is  dependent 
upon  the  predominant  data  trend  of  the  original  data  series,  particularly  in 
the  first  and  last  N  values.  Since  the  low  frequency  trend  in  most  of  the 
lines  of  survey  was  near  linear,  the  data  was  folded  and  inverted.  This 
tends  to  continue  any  linear  trend  present  in  the  beginning  and  end  of  the 
data  and  appeared  to  introduce  the  minimum  error  at  the  beginning  and 
end  of  the  data. 
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APPENDIX  II 


POWER  SPECTRAL  DENSITIES  FOR 
TWO  METHODS  OF  COMPUTATION 
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APPENDIX  III 


POWER  SPECTRAL  DENSITIES 
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Power  Spectral  Density  of  Carswell,  13  R,  filtered  at  500  ft 
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Figure  Id  Power  Spectral  Deusity  of  George*  SOL*  filtered  at  SOO  ft 
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Figure  <  i  Power  Spectral  Density  ol  George,  Center,  filtered  at  500  it 
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Figure  74 


Power  Spectral  Density  of  Glasgow,  Center,  filtered  at  SOD  ft 
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Figure  78  Power  Spectral  Density  of  Longley,  SOR*  filtered  at  S0O  it 
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Figure  T9  Power  Spectral  Density  of  McGuire,  “SL,  filtered  at  500  ft 
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Figure  87 


Power  Spectral  Density  of  Show,  SOU*  filtered  at  800  ft 
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Figure  Si  Power  Spectral  Density  of  Torrejeu,  1 2L,  filtered  at  SOO  ft 
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Figure  SO  Power  Spectral  Deaeits?  o£  Torarejea,  12R.  filtered  at  500  ft 
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Figure  Power  Spectral  Density  of  Travis,  ?§L,  filtered  at  §00  ft 
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Figure  92  Power  Spectral  Density  of  Travis,  Center,  filtered  at  500  ft 
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APPENDIX  IV 


RMS  HISTORIES  FOR  ALL  LINES  OF  DATA 
FILTERED  AT  500  FEET 
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Figure  95  RMS  Histories  for  Carswell  and  Edwards  Filtered  at  500  ft 


Figure  97  RMS  Histories  for  Langley  and  McGuire  Filtered  at  500  ft 
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Figure  96  RMS  Histories  for  Shaw  and  Travis 
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Figure  RMS  Histories  for  Torrejon  and  Castle  Filtered 
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APPENDIX  V 


CUMULATIVE  DISTRIBUTIONS  OF  AMPLITUDES 
FOR  ALL  LINES  OF  DATA  FILTERED  AT  J.00  FEET 
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Figure  104  Cumulative  Distrib/utiea  oi  Amplitude#  tr©r»  Diward* 
Filtered  at  100  it 
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Figure  113  Cumulative  Distribution  of  Amplitudes  from  Torrejon 
Filtered  at  100  ft 
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Figure  114  Cumulative  Distribution  of  Amplitudes  from  Travis 
Filtered  at  100  ft 
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APPENDIX  VI 


RMS  HISTORIES  FOR  ALL  LINES  OF  DATA 
FILTERED  AT  100  FEET 
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Figure  115  RMS  Histories  for  Altus  and  Bakaiar  Filtered  at  100  ft 


m 


SWH 


*  afie«ttituv  SW« 


Ellsworth  and  Glasgow  Filtered  at  100  ft 


Fiftira  l  id  F,M.c  FiV  jries  for  Langley  and  McGuire  Filtered 


Castles 


Figure  121  RMS  Histories  for  George  and  Palmdale  Filtered  at  100  ft 
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